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ABSTRACT 

We present the results of a 10.5 yr, volume limited (28Mpc) search for supernova 
(SN) progenitor stars. In doing so we compile all SNe discovered within this volume 
(132, of which 27% are type la) and determine the relative rates of each sub-type from 
literature studies. The core-collapse SNe break down into 59% II-P and 29% Ib/c, 
with the remainder being lib (5%), IIn(4%) and II-L (3%). There have been 20 II-P 
SNe with high quality optical or near-IR pre-explosion images that allow a meaningful 
search for the progenitor stars. In five cases they are clearly red supergiants, one 
case is unconstrained, two fall on compact coeval star clusters and the other twelve 
have no progenitor detected. We review and update all the available data for the host 
galaxies and SN environments (distance, metallicity and extinction) and determine 
masses and upper mass estimates for these 20 progenitor stars using the STARS stellar 
evolutionary code and a single consistent homogeneous method. A maximum likelihood 
calculation suggests that the minimum stellar mass for a type II-P to form is m m i n = 
8.5~t\ 5 M and the maximum mass for II-P progenitors is m max = 16.5 ± 1.5M0, 
assuming a Salpeter initial mass function holds for the progenitor population (in the 
range T = — 1.35^q7). The minimum mass is consistent with current estimates for the 
upper limit to white dwarf progenitor masses, but the maximum mass does not appear 
consistent with massive star populations in Local Group galaxies. Red supergiants 
in the Local Group have masses up to 25M Q and the minimum mass to produce 
a Wolf-Rayet star in single star evolution (between solar and LMC metallicity) is 
similarly 25-30M Q . The reason we have not detected any high mass red supergiant 
progenitors above 17M0 is unclear, but we estimate that it is statistically significant 
at 2.4(7 confidence. Two simple reasons for this could be that we have systematically 
underestimated the progenitor masses due to dust extinction or that stars between 
17-25M© produce other kinds of SNe which are not II-P. We discuss these possibilities 
and find that neither provides a satisfactory solution. We term this discrepancy the 
"red supergiant problem" and speculate that these stars could have core masses high 
enough to form black holes and SNe which are too faint to have been detected. We 
compare the 56 Ni masses ejected in the SNe to the progenitor mass estimates and 
find that low luminosity SNe with low 56 Ni production are most likely to arise from 
explosions of low mass progenitors near the mass threshold that can produce a core- 
collapse. 

Key words: stars: evolution - stars: supergiants - supernovae: general - supernovae: 
individual - galaxies : stellar content 



1 INTRODUCTION 

Stars which are born with masses above a critical thresh- 
old mass of around 8M0 have long been thought to produce 



supernovae when their cores collapse at a point when nu- 
clear burning no longer provides support against gravity. 
Supernovae were first suggested to be a new class of astro- 
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physical phenomena by iBaade fc Zwickvl fl934) and since 
then detailed study has allowed them to be split into phys- 
ical types : the thermonuclear explosions and core-collapse 
supernovae (CCSNe). The CCSNe form when their cores 
evolve to iron white dwarfs and detailed stellar evolution- 
ary models pred ict a minimum mass for this to occur of 
between 7- 1 2MrrjHeger et al]l2003l; lEldridge fc Toudl2004bl ; 



Sicss] 120071 ; iPoelarends et al 



2008) . The CCSNe form a 



diverse group in terms of their spectral and photometric 
properties and the classification scheme has arisen primar- 
ily based on the appearance of their optical spect ra, but also 
supp lemented with their photometric behaviour (|Filippenkd 
1 19971 ). There are the H-rich type II SNe, which are sub- 
classified into II-P (plateau lightcurves) , II-L (linear decline 
lightcurves) , Iln (narrow emission lines), and some peculiar 
events, generically labelled II-pec. The H-deficient SNe are 
split into lb and Ic depending on whether He is visible and a 
hybrid class of lib (type II events which metamorphose into 
lb SNe) has also been uncovered. The evolutionary stage of 
the progenitor star (i.e. its position in the HR diagram and 
chemical composition of its atmosphere) very likely dictates 
what type of CCSN is produced. 

Determining what types of stars produce which types 
of SNe and what mass range can support a CCSN is a ma- 
jor goal in modern studies of these explosions. The first at- 
tempts relied primarily on linking the results of stellar evo- 
lutionary models to the spectral (and photometr ic) evolu- 
tion o f the SNe (e.g . lArnettlll980l;lchevalierlll976l ) More re- 
centlv lHamuvl (|2003 ) and lNadvozhinl |2003h have studied the 
lightcurves and ejecta velocities of II-P SNe to estimate the 
mass of the material ejected. These studies tend to favour 
quite large masses for progenitor stars, with Hamuy sug- 
gesting ranges of 10-50Mq and Nadyozhin 1O-3OM0. How- 
ever when SN 1987A exploded a new opportunity arose. 
The SN was in a galaxy close enough that its progenitor 
star could be easily identified. A blue supergiant star of 
around 15-2OM0 was i dentified and it is clear that this ob- 
ject no longer exists (iGilmozzi et al.l 1 19871 ; IWalborn et al.l 
1 19891 ; |Podsiadlowskilll992l ) The fact that it was a compact 
blue supergiant was a key factor in enabling the community 
to understand the event as a whole. A progenitor was also 
detected for the next closest explosion (SN 1993J in M81) 
and the binary nature of the progenitor helped understand 
the physical reason behind the H-b type given to the event 
jAldering et all 1 1994 iPodsiaklowski et alj|l993h . 

Studies of the environments of CCSNe after discov- 
ery have been ongoing for many years with authors look- 
ing for correlations between the ages of starform ation re- 
ions and t he type of SNe t hat o ccur. For example Ivan Dvkl 
19921 1 and Ivan Dvk et all (|l99d l suggested that there was 
no clear difference in the spatial distributions of type Ib/c 
and type II SNe compared to giant Hn regions in their 
host galaxies. They concluded that they hence arose from 
parent populations of sim ilar mass. Howeve r recen tly both 
iJames fc Andersonl (|2006l ) and iKellv et al.l (|2008l ) suggest 
that the stripped type Ic SNe are more likely to follow re- 
gions of either high surface brightness or high Ha emis- 
sion than type II SNe. A larger sample of nearby core- 
collapse events and Ha c orrelations has been compiled by 
l|Anderson fc Ja mes 2008) indicating that there is a progres- 
sive trend for Ib/c SNe to be associated with Ha emission 
regions, in the sense that Ic show the closest association with 



galactic Ha emission, then comes the lb and then type II. 
While these efforts are very valuable to discern differences 
in progenitor channel, it is difficult to assign definitive mass 
ranges to the progenitor systems from spatial correlations 
alone. 

A much more direct way to determine the type of star 
that exploded is to search directly for progenitors in im- 
ages of the host galaxies taken before explosion. The ease of 
access to large telescope data archives makes this search fea- 
sible for nearby events. There are now a number of groups 
around the world that are competitively searching for pro- 
genitor stars in such archive images. Particularly with the 
Hubble Space Telescope (HST) it has become possible to re- 
solve massive stellar populations out to at least 20Mpc. In 
this case the information on each progenitor is much more 
detailed and quantitative than can be achieved with the 
unresolved environment studies, but there are fewer events 
which allow such a study. 

Early work wit h HST concentrated on looking at the en - 
vironments of SNe l|Barth et al.lll99rj ; IVan Dvk et al.lll999l '). 
But it has now become possible to search directly for progen- 
itor stars and carry out concerted campaigns on the nearest 
events. Three of the first type II SNe to have excellent HST 
and ground-based pre-explo sion images were the II-P SNe 
1999gi, 1999em and 2 001du (|Smartt et al.ll200ll . l2002l . l2003l ; 
I Van Dvk et~ai1l2003bf ). In all there was no det ection of a pro- 
genito r, bu t meaningful limits w ere derived. iLeonard et al.l 
( 200$ and ILeonard et all (|2002h subsequently showed the 
importance of having reliable distances to the host galaxies 
of the SN progenitors in order to constrain th e upper mass 
limit s . Efforts to find proge nitors continued (|Smartt et al.l 
120031 ; IVan Dvk et all 1200331 ) until the first confirmation of 
a red supergiant progenitor of a type II-P explosion in 
predi scovery HST and Gemini-North images (|Smartt et alj 
|200J). Unambiguous detections in HST images require the 
SN to be located on the pre-explosion images with accura- 
cies of around 10 milliarcseconds, which requires follow-up 
images of the SN to be taken at HST or corrected adap- 
tive optics ground-based resolution ( Maund fc Smartt]|2005l ; 
iGal-Yam et al.l 120051 ; ICrockett et all |200ST ) The next clear 
and unambiguous detection of a progen itor star was also a 
red supergiant in NGC5194 (SN2005cs ; iMaund et ai1l2005l ; 
iLi et al.ll2006h . And the recent discovery of SN2008bk in 
NGC7793 (3.9 Mpc) has produced detec tions in IJHK of a 
red progenitor star l|Mattila et al.ll2008r ). The near-IR spec- 
tral energy (SED) distribution of the SN2008bk progenitor 
is the best sampled SED yet of any red supergiant progen- 
itor and matches a late M4I spectral type with moderate 
extinction of Av — 1 and an initial mass around 8-9Mq. 
There have been claims of the dete ctions of other s such as 
2004A (|Hendrv et al"]|2006l) , 2004et (|Li et al.ll2005l ). 2006my 
and 2006ov (ILTet al. 2007J) and we review these in this pa- 
per. An additional method that has much to offer this field 
is locating SNe directly coincident with compact, coeval star 
clusters. If the age of the cluster can be determined then a 
main-sequence turn-off age, and turn-off mass can lead di- 
rectly to an estimate of the progenitor star mass (e.g. 2 004dj 
in NCG2403 as studied by iMafz- Apellaniz et al"1l2004l '). The 
detection and characterisation of progenitor stars has the 
potential to directly link the type of star to the SN explo- 
sion characteristics, (e.g. the amount of 56 Ni synthesised in 
explosive O and Si burning and the total energy of the ex- 
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plosion) and also to set quantitative limits on progenitor 
mass ranges. The latter is of great interest to compare to 
the highest mass white dwarf progenitors and to the stars 
that form neutron stars and black holes after core-collapse. 

Work in this field has progressed substantially in the 
last 8 years an d compilation s of pr o genitor properties hav e 
been made bylSmartt et al.1 ll2003h: leal-Yam et all (|2007h : 
iLi et all i|2007yKochanek et al.1 (|2008T ) . However these com- 
pilations are somewhat ad-hoc, incomplete and potentially 
biased as they do not define the selection criterion for inclu- 
sion rigorously. In addition the results are based on differ- 
ent methods for estimating the progenitor masses and upper 
limits (in terms of measurement and the theoretical models 
employed) . The goal of this series of two papers is to define 
the selection criteria for inclusion (a volume and time lim- 
ited survey) and to determine the physical parameters of the 
progenitor stars (luminosities and masses, or limits thereon) 
in a homogeneous and consistent way. Only then is it pos- 
sible to reliably estimate the population parameters. This 
paper specifically deals with the type II SNe and all but two 
of the final sample of twenty have been confirmed as type 
II-P. A companion paper will discuss the stripped SNe of 
types lib, lb and Ic which are drawn from the same volume 
and time limited survey (Crockett et al. 2008, in prep). 

In addition to the supernovae with pre-explosion images 
we discuss in the paper there are three SNe with progeni- 
tor detections that fall outside either our time or volume 
definition. Those are SN 1987A, SN 1993J and SN 2005gl. 
Discussions of the first two are we l l documented i n the 
literature (e.g IWalborn et al.1 Il989l; lAldering et al.1 1 1994 



I Van Dvk et~aLll2002l ; iMaund et al.ll2004r i; 1987A will be dis- 
cussed later in this paper and the implications of 1993J will 
be discussed in the 2nd paper in this series. SN 2005gl is 
a type Iln in NGC266 at approxi mately 66 Mpc, hence al- 
though a progenitor is detected bv lGal-Yam et aL (12007) it 
does not fall within our distance limit. iGal- Yam et al. I (120071 ) 
suggest that this was a very massive star with Mv — —10.3 
and very likely a luminous blue variable. Although with only 
a detection in a single filter the blue colour is not confirmed, 
and there is not enough data to determine if the source was 
indeed variable. This is evidence that very massive stars do 
explode as bright SNe and is a point we will return to in the 
discussion. 

This paper starts with defining the sample from which 
the targets with high quality pre-explosion images are 
drawn. A consistent and homogeneous analysis method is 
then defined and justified. We then review previous detec- 
tions (and add some new data) to build the data required 
for analysis and then statistically analyse the results. We 
follow this with an extensive discussion. 



2 THE SAMPLE OF LOCAL UNIVERSE 
SUPERNOVAE 

The observational data for this paper are compiled from 
many sources in the recent literature but the sample selec- 
tion requires some justification and explanation if the later 
comparisons and discussions of physical parameters are to 
be meaningful. We have selected SNe for inclusion based on 
the following selection criteria, and we justify the choice of 
criteria where appropriate. 



2.1 Definition and selection 

We consider all core collapse SNe discovered in the ten year 
period between 1998 January 1 and 2008 June 30. The ear- 
lier date was chosen as the sensible start point for the con- 
certed efforts to find SNe in archive pre-discovery images due 
to the fact tha t the local SN discovery rate had reached a sig- 
nificant level (Ivan den Bergh et al.1 120051 ). and we estimate 
the amount of imaging of nearby galaxies in the Hubble Space 
Telescope (HST) archive had become rich enough that coin- 
cidences were likely to occur. This was also the effective date 
of the start of concerted efforts to search for SN progenitor s. 
Since the late 1 990's our group Ipmartt et alj|200ll. |2002|) 
and others (e.g. lVan Dvk et al.ll2003al ; iGal-Yam et alll2005l ) 
have been systematically searching for pre-explosion images 
of core-collapse SNe in nearby galaxies. We have further re- 
stricted our sample to galaxies with recessional velocities 
less than 2000 kms -1 , effectively restricting us to a volume 
limited sample. The recessional velocities, corrected for the 
infall of the Local Group towards Virgo cluster, of all nearby 
galaxies hosting SN e were taken from the HyperLEDALj 
l|Paturel et al.l [20031 ) database. We emphasise that we have 
used the corrected velocities to apply the selection criterion 
of 2000 km s~ 1 and assuming Ho = 72km s~ 1 Mpc - 1 , this lo- 
cal volume has a limit of 28 Mpc. 

During this period, and in this volume, there have been 
138 SNe candidates discovered and all are listed in the Asi- 
agdl CfA0 

and Sternberg Astronomical Institute |_| cata- 
logues. Of course these are simply compilations of discover- 
ies reported in the International Astronomical Union (IAU) 
Circulars and these catalogues normally list the discovery 
magnitudes and types reported in the first IAU announce- 
ment. It often happens that the SN classification is revised 
or refined with subsequent higher quality spectra, or longer 
monitoring. Both can reveal peculiarities and transforma- 
tions, or simply give a more secure classification. In partic- 
ular the sub-classification of II-P can be added when signifi- 
cant light curve information is gathered. Hence we have care- 
fully checked the classification of each event, and have gone 
further and classified those supernovae listed as "II" into 
the subtypes II-P and II-L or Iln where possible. This was 
done using the following criteria in order. Firstly the refereed 
literature was searched and a classification taken from pub- 
lished photometric and spectroscopic results. Secondly un- 
published, professional spectra and lightcurves where taken 
from reliable so urces such as those of the Carnegie Su- 
pernova Project iHamu^etj^ [200(^) and the Asiago data 
archive l|Turattoll2000l ). Thirdly amateur lightcurves avail- 
able on the web were checked and, if possible, a II-P classifi- 
cation was made if a clear and unambiguous plateau lasting 
longer than 30 days was recorded. The vast majority of II-P 
SNe have plateau phases lasting significantl y longer than 30 
days, in fact most are around 90-110 days 1 Pastorelioll2003l ; 
iHamuvl 120031 ; iPastorello et al.l 120041 ') and there is no clear 
evidence that plateaus of shorter duration are particularly 
common. However to observe such a long plateau necessarily 
means the SN must have been discovered close to explosion. 



1 http://leda.univ-lyonl.fr 

2 http:/ /web. pd.astro.it/supern/snean.txt 

3 http:/ /cfa-www. harvard.edu/cfa/ps/lists/Supernovae. html 

4 http: / /www. sai.msu.su/sn/sncat/ 
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This is often not the case, and we have chosen to take 30 
days simply as an indicator that an extended plateau phase 
is evident. In all cases of our sub-classifications of type II 
SNe, we believe the designations not to be controversial or 
ambiguous and for only 9 events classed as type II were we 
unable to assign a subtype. In these cases the supernovae 
were generally discovered late in the nebular phase. It is of- 
ten difficult to distinguish lb and Ic SNe from single spectra 
taken at an unknown epoch, and indeed there are 6 events for 
which authors have listed Ib/c classifications and we cannot 
improve on these classifications. All of the individual SNe 
are listed in Table lATI There are only two SNe which have 
not had a classification spectrum reported in the literature, 
1998cf and 1999gs and these are ignored in the following 
frequency comparison. 

We realise that the classifications into the standard 
bins is somewhat simplistic. In particular there are some 
SNe that show evidence of interaction with the circumstel- 
lar medium, which result in narrow lines (usually of H or 
He) superimposed on the spectrum. When narrow lines of 
H dominate the spectrum then the Iln designation is of- 
ten used by the community, but some II-P and Ibc SNe 
do show evidence of this behaviour at a weaker level. The 
most striking example recently is that of the SN2006jc-like 
objects and the two examples in our sample are 2006,jc 
and 2002ao (see IPastorello et all 120071 ; iFolev et al.l 120071 ; 
iPastorello et al.ll2008l ). These show broad lined spectra re- 
sembling type Ic SNe (in that they do not exhibit H or 
He in the high velocity ejecta), but have strong and nar- 
row He emission lines and weak Ha emission. This has led 
IPastorello et al.l (|2008h to term this class of objects Ibn. 
However rather than introducing this small and very spe- 
cific class of events, we will class them as Ic SNe as this is a 
fair description of the underlying spectrum of the SN ejecta. 
Giving them a simple lb label could be argued as being mis- 
leading in that they do not exhibit the broad, He absorption 
typical of this class of H-deficient events. The label Ibn is 
certainly a valid type for them but in this paper it is too 
specific to be a useful addition to the compiled subtypes. 
We will discuss these objects further in the second paper in 
this series which concerns the stripped events (Crockett et 
al. in prep.). 

2.2 The relative frequencies of core-collapse SNe 

In TablefTJwe list the relative frequency of each subtype oc- 
curring in our sample. This is a volume limited relative fre- 
quency rate of SN types. There may well have been undiscov- 
ered local SNe in this period e.g. dust extinguished events, 
or events which exploded in solar conjunction which were 
missed at late times. The distance limit imposed (/x=32.3) 
and the range in the ab solute magnitudes of each subtype 
l| Richardson et all |2002| ) would initially suggest that it is 
unlikely that there is a serious bias in the relative num- 
ber of the different subtypes. However there are two argu- 
ments which can be put forward against this. The first is 
if there is a substantial number of intrinsically faint SNe 
that may have gone undetected if they have typical mag- 
nitudes below ab out —13. The nature of the faint opti cal 
transient in M85 (|Kulkarni et al.ll2007l ; Ibfek et al.ll2008l ) is 
curr ently debated, and co uld conceivably be a core-collapse 
SN l|Pastorello et al.ll2007l . see Sect. 2.4). If a large number of 



intrinsically very faint SNe are evading discovery by current 
pointed surveys it could lead to a dramatic change in our 
understanding of the link between progenitor star and SN 
explosion. Secondly there is an issue with the lack of type la 
SNe discovered recently within lOMpc. In the 10.5 yr con- 
sidered there are about 13 core-collapse that have been dis- 
covered within 10 Mpc (the exact number depends on some 
individual galaxy distance estimates and how strictly one 
enf orces the distance limit; for a more in depth discussion 
see lKistler et al.ll2008l ). But there have been no type la SNe 
discovered and with a relative frequency of 27%, one might 
have expected to have seen 3-4 SNe. The Poisson probabil- 
ity of this being a statistical fluctuation is not zero, but is 
small (2-5%) and one could invoke the argument that there 
are more core-collapse SNe (perhaps of the fainter type II) 
beyond 10 Mpc which are being missed and hence the rela- 
tive rate of CCSNe/Ia is i ntrinsically much high er than we 
curr ently believe (also see Thomps on et alj|2008f ). 

Ivan den Bergh et al. I (120051 s ) have presented a homoge- 
neous sample of 604 recent SNe discovered (or recovered) 
by the Lick Observat ory Supernova Search (LOSS) with 
the KAIT telescope i|Filippenko et all I2OO1T ). The galaxy 
search sample spans a much larger volume (cz < 10 000 
kms -1 ) than we are considering and a significant major- 
ity of our sample in the overlapping time frames (~85 per 
cent of the core collapse events between 1999-2004) are 
listed in the van den Bergh et al. summary of the LOSS 
survey. The ones which are not are predominantly more 
southern than 5 = —30°. Hence our sample is very simi- 
lar to that which would be obtained if one selected a dis- 
tance and time li mited sample from the disc overed and re- 
covered events of Ivan den Bergh et all (|2005l ). The relative 
number of type la SNe in the full LOSS catalogue is sig- 
nificantly higher than within our smaller volume (44 per 
cent compared to 27 per cent that we find here). This is 
very likely due to type II-P SNe going undetected at the 
largest distances. At cz ~ 10000 kms -1 , with a moder- 
ate amount of fo reground reddening at le ast half of the II- 
P distribution of iRichardson et al.l (|2002l ) would be missed 
at the limiting magnitude of R ~19 of the KAIT survey. 
The type Ib/c SNe appear slightly more abundant in our 
loc al sample (29 per cen t of al l core-collapse) compared to 
the Ivan den Bergh et al.l (|2005h frequency (25 per cent), al- 
though the difference is not significantly greater than the 
expected Poisson scatter. 

The 9 SNe which were classed as type II (and could not 
be further sub-classified) can be split proportionally over 
the type II subtypes, which assumes that that there was no 
particular bias underlying their poor observational coverage. 
As they were all discovered late in the nebular phase, this is 
likely to hold. We put 8 in the II-P bin, and split the other 
1 equally between the Iln and ILL, as they have equal num- 
bers of confirmed types; hence the fraction which appears 
in table [T] In a similar manner the 6 Ib/c SNe are split 
proportionately into the lb and Ic bins based on the mea- 
sured ratio of Ib:Ic = 7:14 (which comes from those events 
where a lb or Ic classification seems secure). We did not in- 
clude the 2 unclassified SNe in any of the rate estimates. 
Seven of the events originally announced as SNe in Table 
I All have been shown to actually be outbursts of luminous 
blue variables (LBVs) similar to those seen historically in 
Local Group LBVs such as 77-Carinae and P-Cygni. These 
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Table 1. The relative frequency of SNe types discovered between 
1998-2008 (10.5 yrs) in galaxies with recessional velocities less 
than 2000 kms -1 , and type taken from Table lAll The relative 
frequency of all types and the relative frequency of only core- 
collapse SNe are listed separately. 







Relative 


Core-Collapse only 


Type 


No. 


/ per cent 


/ per cent 


II-P 


54 


39.1 


58.7 


II-L 


2.5 


1.8 


2.7 


Iln 


3.5 


2.5 


3.8 


lib 


5 


3.6 


5.4 


lb 


9 


6.5 


9.8 


Ic 


18 


13.0 


19.6 


la 


3V 


26.8 




LBVs 


7 


5.1 




Unclassified 


2 


1.1 




Total 


138 


100 


100 


Total CCSNe 


92 


66 


100 



are 1999bw jFilippenko et al.lll999l). 2 000ch (|Wagner et all 
l2004h , 2001ac jMatheson fc Calkins} 1200 ll). 2002kg or 
NGC2403-V37 JWeis fc Bomansl 12003; flVIaund et alj 120061 ; 
IVan Dvk etHu1|2006l) . 2003gm dMaund et al.ll200rj) 2006f p 
l|Blondin et~aTT l2006h . 2007sv (jHarutvunvan et all l2007h . 



Hence we remove these from the rates of CCSNe since they 
are not true SN explosions. Table [1] lists the relative fre- 
quencies of core-collapse SNe. It is clear that the types II-L 
and Iln are intrinsically quite rare and the majority of core- 
collapse events are SNe II- P. Such a breakdown of subtypes 
has b een suggested before (|Li et al.ll2007i ; ICappellaro et al.l 
1 19991 ) although this is the first time quantitative volume lim- 
ited statistics have been compiled and presented. The pre- 
liminary analysis bv lLi et"a!1 (|2007l ) of 68 LOSS only discov- 
ered events (within 30Mpc) in 9 years suggests 68:26:2:4 per 
cent breakdown between II:Ib/c:IIb:IIn. Perhaps the ratio 
of most interest is the Ibc/II ratio which has been used by 
previous stu dies to try to plac e constraints on progenitor 
populations ( Prieto et al.1 120081 : iPrantzos fc Boissierl 120031: 
lEldridgdl2007f r These three studies have estimated the ratio 
as a function of metallicity finding that, at approximately 
solar metallicity (Zq), the ratio is Nn, c /Nu ~ 0.4±0.1 (Pois- 
sonian uncertainty) and this goes down to around 0.1 at 
O.3Z0 (although with fairly small numbers in each metallic- 
ity bin). The ratios at approximately Z© are fairly similar 
to what we find (Nj h / C /Nu = 0.45 ± 0.13) and as discussed 
below in Sect(4]and Sect[5]our SN population is likely drawn 
from metallicities in the range O.5-1.OZ0 due to the fact that 
nearby, high starformation rate galaxies are those that are 
most frequently monitored for SNe. A full analysis of the 
chemical composition of the sites of the SNe in this volume 
limited sample would be desirable. We will discuss the SN 
rates more in Sect. 18.21 

The absolute rates of the different types in the stan- 
dard supernova units (ISNu = 1 SN(100yr)- 1 (10 10 L|)- 1 ) 
is much more difficult to assess given that detailed knowl- 
edge of the sampling fr equency for each galaxy and search 



strategy is required (see Cappcllar o et al 



team will address this l|Leaman et al.l 2004 ) and will provide 
the best estimate of the local rates so far. While we cannot 



19991). The LOSS 



derive the SN rate in standard units, the numbers in Ta- 
ble[T] serve as a good guide to the relative numbers of SNe 
expected in future surveys, when used in conjunction with 
absolute magnitude distributions (assuming the local galaxy 
population is cosmically representative). If bias factors af- 
fecting the relative rate of discovery of the different types 
are minor, then these rates are a direct consequence of the 
initial mass function combined with stellar evolution which 
is dependent on mass, metallicity, duplicity and initial ro- 
tation rate. The numbers of SNe we have tabulated give a 
lower limit on the number of types per Gpc _3 yr~ 1 , which 
is at least a useful comparison to higher redshift estimates 
and also when considering the rates of Gamma-Ray Bursts 
(GRBs) and X-ray Flashes (XRFs). The volume enclosed by 
the 28Mpc limit is 9.2 x 10~ 5 Gpc 3 ; hence the local rate of 
CCSN explosions is likely to be > 9.6 x 10 4 Gpc^yr -1 and 
the local Ibc SN rate is ^ 2.7 x 10 4 Gpc~ 3 yr~\ The latter 
is in reasonable agreement with the ~ 2 x 10 4 Gpc _3 yr _1 
put forward by iGuetta fc Delia Valid (2007), based on the 
ICappellaro et al.1 (|l999F T rates and local galaxy luminosity 
functions. 

One further point to bear in mind when considering 
the relative rates is the likely number of local SNe which 
are not discovered because they are in faint hosts which 
are not monitored. The Sloan Digital Sky Survey Data Re- 
lease 5 (5713 sq degrees) contains about 2200 faint galax- 
ies with M g —17 and with recessional velocities less 
than 2000 kms -1 . Hence over the full 40,000 square de- 
grees of sky, there are about 15,000 of these faint hosts 
within about 28 Mpc. These are generally not monitored 
by the LOSS and the amateur efforts, who typically tar- 
get the most luminou s 10,00 galaxies within about 100- 
140 Mpc. lYoung et al.l ((2008) estimate that such galaxies 
(with metallicities corresponding roughly to oxygen abun- 
dance < 8.4 dex) would contribute about 5-20% of the total 
star formation locally. Hence at least this fraction of core- 
collapse SNe are missing from the local samples and all of 
them are within faint hosts and low metallicity progenitors. 
This could mean t hat very bright events like SNe 2005ap , 
2006tf and 2008es (lOuimbv et"ai1l2007l : ISmith et al1l2008al : 
iGezari et all 120091 : [Miller et al.l 120091 ) found in blank-field 
searches (and faint hosts) could be missed. Although the 
true rate of suc h events appears t o be quite small and likely 
less than ~1% (Miller et al. 2009). Discovery of these events 
locally should be possible with future all-sky surveys such 
as Pan-STARRS and LSST (for an estimate of rates see 
lYoung et al"1l2008r i. 

Of the 100 core-collapse SNe and LBV classified out- 
bursts, the host galaxies of 46 of them were imaged by HST 
with either the Wide-Field-Planetary-Camera 2 (WFPC2) 
or Advanced Camera for Surveys (ACS) before explosion. 
However given the small field of view of both of these cam- 
eras (2.6' and 3' respectively) the site of the SNe did not al- 
ways fall on the field-of-view (FOV) of the cameras. Of these 
46, only 26 had images of the SN site in the camera FOV, 
an overall hit rate of 26 per cent. A column is included in 
Table lATI which specifies whether or not the galaxy was ob- 
served by HST before explosion and if the SN falls on one of 
the camera FOVs. A further 6 have had high quality ground- 
based images of the SN site taken before explosion and they 
are also included in this compilation. The observational sam- 
ple for this paper, and its companion studying the stripped 
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events (Crockett et al. 2008, in prep.) is thus all of the core- 
collapse SNe which fulfill the above criteria and have good 
quality pre-explosion imagery. We have confirmed that there 
are no other SNe in Table lXTI with HST pre-discovery im- 
ages. We cannot make the same definitive statement about 
ground-based images given the amount of inhomogeneous 
imaging data around. But our manual searching of all well 
maintained large telescope archives suggests it is highly un- 
likely that further high quality images of any of these events 
will surface i.e. images with sub-arcsec resolution with the 
depth to detect a large fraction of the galaxy's massive stel- 
lar population. As such we have a well defined sample in 
terms of distance and time. The rest of this paper focuses 
on the progenitor properties of the 20 type II SNe listed in 
Tabled of which 18 are confirmed II-P and two are of un- 
certain subtype (1999an and 2003ie). The other twelve su- 
pernovae which are likely to have had stripped progenitors: 
2000ds (lb), 2000ew (Ic), 2001B (lb), 2001ci (Ic), 2002ap 
(Ic), 2003jg (Ib/c), 2004gt (Ib/c), 2005ae (lib), 2005V (Ibc), 
2005cz (lb), 2007gr (Ic), 2008ax (lib), will be discussed in a 
companion analysis paper (Crockett et al. in prep.). 

2.3 The relative frequencies of SN1987A-like 

events 

In this volume limited sample, there is only one SN which 
has been con clusively shown to be s imilar to SN1987A, that 
is SN1998A (|Pastorello et all 120051 ). SN1987A had a pecu- 
liar lightcurve and distinctly strong Ball lines (probably a 
temperature effect) and an asymmetric Ha profile during 
its first ~40 days of evolution and the community would 
have been very unlikely to miss such events as they would 
have created great interest. There is one other relatively 
nearby ev ent that has a SN1987A-like appearance, which 
is 2000cb i|Hamuvll206"iT ) in IC1158. This one however has 
a V v ir = 2017 kms -1 , which puts it jus t beyond our selec- 
tion cr iteria. We shall see in Sect l5.9l that lHarutvunvan et al.l 
(2008) suggest that a single spectrum of SN2003ie shows 
similarities to 1987A but it is not well studied enough to 
be definitive. Hence even if we would include 2000cb and 
2003ie in our sample we can certainly say that 1987A events 
are intrinsically rare and probably less than around 3% of 
all core-collapse events. 

2.4 SN2008S and the M85 and NGC300 optical 
transients 

Recently three optical transients have been reported whose 
nature is still ambiguous and in tensely debated . The opti- 
cal transient in M85 reported bv lKulkarni et al.l (|2007h was 
suggested by the authors to be a "luminous red nova" which 
most likely arose from a stellar merger. H owever this view 
was challenged by iPastorello et al.l 1)20071) who suggested a 
CCSN origin could not be ruled out. Since then two other 
optical transients of similar absolute magnitude have been 
discovered. One has been ter med a supernova (S N2008S; 
IStanishev et al.l 120081 ) although ISmith et all (|2008l ) suggest 
it could be a supernova imposter and the outburst of a mod- 
erately mass ive star rather than a core-collapse. T he other, 
in NGC300 <|Bond et al.l 120091 ; IBerger et af]|2009l ). has not 
yet received an official supernova desig nation, hence we re - 
fer to it as NGC300 OT2008-1 (as in IBerger et~aH |2009| ) 



SN2008S has already been subject to a study of its pre- 
explosion environment and a detection of a source in Spitzer 
mid-IR images has been suggested to be a dus t enshrouded 
red supergiant which is visually obscured |Prieto et al.l 
2008). A similar dust dominated object has been found to 
be coincident with the optical transient NGC300 OT2008-1 
llThompson et al.ll2008l : iBond et al.ll2009l ). iThompson et alj 
(2008) suggested that all three could be the similar explosion 
of massive stars embedded in optically thick dust shells. The 
early studies of the evolution of SN2008S, NGC300 OT2008- 
1 and their comparisons with M85OT2006-1 and other 
erupting systems have so far not favoured a core-collapse su- 
pernova explanation for the physical source of the out burst 
l|Smith et al.ll2008l ; lBond et al.ll2009l ; IBerger et al.ll2009l ) The 
transients lack broad lines from high velocity ejecta; their 
spectra are very slowly evolving and dominated by narrow 
H-emission. Strangely they also don't appear to be similar 
to the V838 Mon variable system or M31 luminous red vari- 
able as initi a lly su ggested by both lKulkarni et al.l (|2007l ) and 
IBond et"aL ||2008 | ). Ba s ed on the mid-IR progenitor detec- 
tions. lThompson et al. I l|2008l) argue that the precursors may 
have been going through a short evolutionary phase which 
ends in a weak, electron-capture supernova. A full multi- 
wavelength study of the evolution of SN2008S from early to 
late times, and comparisons with the other two suggest there 
is some evidence for the supernova explanation (Botticella 
et al., in prep.; Kotak et al. in prep). All these studies reveal 
that the three objects are incredibly similar in their proper- 
ties. As their nature is ambiguous and currently debatable, 
we will not consider them further in this paper. It is certain, 
however, that they are not normally type II-P SNe. 



3 THE STELLAR EVOLUTIONARY MODELS 

As discussed above, observational and theoretical studies 
both now strongly suggest that the progenitors of type 
II-P are typically red supergiants. To estimate an initial 
mass for observed red supergiant progenitors we require 
stellar models to obtain a theoretical initial mass to fi- 
nal luminosity relation, as shown in Figure [1] The stellar 
models we use were produced with the Cambrid ge stellar 
evolut ion code, STARS, originally d eveloped b y Eggletonl 
Jl97ll) and updated mos t recently by iPols et all l| 19951 ) and 
lEldridge fc Toutl (2004a). Further details can be found at 
the code's web pageqj. The models are available from the 
same location for download witho ut restriction. The models 
are the same as those described in lEldridge fe Toutl (|2004al ) 
but here we use every integer initial mass from 5 to 40Mq 
and integer steps of 5-IOMq above. 

As will be discussed in Sect .[4] we can estimate the 
metallicity of the exploding star from the nebular abun- 
dances in the disks of the host galaxies, hence we have cal- 
culated stellar evolutionary models for three metallicities; 
solar, LMC and SMC where we assume mass fractions of 
Z=0.02, 0.008 and 0.004 respectively. All the models employ 
our standard mass-loss p rescription for hydr ogen-rich stars 
l|Eldridge fc Toutll2004bl ): we use the rates of Ide Jager et all 
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(1988) except for OB star s, for which we use the theoretical 
rates of IVink et all l|200ll ). 

In Figure [T] we plot the range of luminosity for a star 
from the end of core helium burning to the model end point 
at the beginning of core neon burning (for a solar metallic- 
ity model). The beginning of core neon burning is only a 
few years before core-collapse and this point is likely to be 
an accurate estimate of the pre-SN luminosity. The estimate 
of final mass from the observational limits will depend on 
uncertainties in these stellar models, which is a systematic 
that is difficult to constrain. To allow for this we assume that 
the range of reasonable luminosities for progenitor stars is 
somewhere between the end of core helium burning (dashed 
line in Fig. la) to the model end point at the beginning of 
core neon burning (solid line in Fig. la). For the lower mass 
stars that undergo the process known as second dredge-up 
to become AGB stars, we also consider the luminosity before 
second dredge-up occurs. After second dredge-up the mod- 
els have much higher bolometric luminosities but their ob- 
servable characteristi cs are quite different to red supergiants 
(|Eldridge et al.ll2007l ). We have previously shown that in the 
case of SN2005cs the progenitor could not have been a super- 
AGB star in the 5-8M0 range. Hence we assume throughout 
this paper that such stars are not the progenitors of the SNe 
discussed, and their positions in Fig.[T] are s hown for com- 
plete ness. A full discussion of this is given in (|Eldridge et al.l 
120071 ). In Fig.QJ) we show the final luminosity ranges for the 
three metallicities, and the most appropriate metallicity for 
each SNe is used when the initial masses are calculated in 
Sect.0 

In Sect. 5, when we estimate a progenitor initial mass 
from Figure [1] we will assume that the models are reliable 
enough to predict correctly that SNe will undergo core- 
collapse after helium burning, so we use the full range of 
luminosities between the start of helium burning and the 
model endpoints in a conservative way. The mass estimate 
should be reasonable for all cases where the progenitor is a 
red supergiant. 



4 METALLICITIES OF THE PROGENITOR 
STARS 

As the stellar evolutionary tracks do differ slightly, to make 
sure that there are no underlying systematics in our anal- 
ysis we require an estimate of the initial metallicity of 
the exploding star. There is good observational evidence 
now to show that mass-loss from massive stars is metal- 
licity dependent, and that the the lifetimes of stars in 
various phase s as they evolve dep end on metallicity (e.g. 
iMassevI 120031 ; iMokiem et~al1 l2007t ). Models predict that 
mass-loss and metallicity are d r iving forces behind stel- 
lar evolution drleger et all 120031 ; iMevnet fc Maederl 120031 ; 
lEldridge fc Toutll2004bl ). 

The most reliable determination of the metallicity of 
the progenitor star would be a measurement of the inter- 
stellar medium abundance in the galactic disk at the posi- 
tion of the event. In some cases the SNe have exploded in, 
or very close to, a previous catalogued Hn region which 
has published spectroscopy and emission line fluxes. For 
these events we use these fluxes and calculate the nebu- 
lar abundance of oxygen, using the strong-line method of 
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Figure 1. (a): The initial mass compared with the final lumi- 
nosity of the stellar models for Z = 0.02. Each mass has the 
luminosity range corresponding to the end of He-burning and the 
end of the model, just before core-collapse (these are the thin grey 
vertical lines). From a limit of luminosity an upper limit to the 
initial mass can be determined. The solid line is the luminosity 
of the model end-point, the dashed line the luminosity at the end 
of core helium burning and the dash-dotted line is the luminosity 
after second dredge-up when the lower mass stars become AGB 
stars, (b): The same as (a) but with three metallicities shown for 
comparison, and with the vertical joining bars omitted for clarity. 
This illustrates that the choice of metallicity for the tracks is not 
critical, but we do use the most appropriate track to remove any 
systematic error. 



the rati o of the [OH] AA 3727 plus [O III] AA 4959,5007 to 
H/3 (see lBresolin1l2006l . for a discussion). Recently there has 
been much debate in the literature over which calibration 
to use to de termine the ne bular oxygen abundances from 
th is method. | Bresolinl |2006| ) has shown that the calibration 
of iPilvugin fc ThuanT l|2005h best matches the abundances 
determined in nearby galaxies where it is possible to mea- 
sure the strength of the electron temperature sensitive lines 
and hence determine a simple empirical calibrations for the 
strong-line m ethod. The result is that the empirical deter- 
minations of Bresolin et alj d2004h and the ca li bratio ns of 
IPilvugin fc Thuanl [j20ol Taiid IPettini fc Pagell (|2004l ) give 
significantly lower abundan ces (by a factor 0.3-0.4 d ex) than 
photoionisation models (of iKewlev fc Dopital |2002t for ex- 
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ample). iTrundle et all (2002) have shown that photospheric 
abundances of massive stars in M31 (B-type supergiants) 
are in much better agreement with th e "lower" metallic- 
ity sca les of the P method employed by IPilvugin fc Thuanl 
(2005), and the empirical determinations form auroral lines. 
The recent downward s revision of the solar oxyg en abun- 
dance to 8.66 ± 0.05 (|Asplund et all 12004 120051 ) and the 
agreement with all the other estimators of the Milky Way's 
ISM abundance at the solar radius (B-stars, young F& G 
stars, Hn regions, diffuse ISM) wou ld also suggest that our 
adopted, lower, scale is appropriate dSimon-Di'az et al . 2006; 
iDaflon et alll2004l : ISofia fc Meverll200ll). Hence in this p aper 
we will favour the calibr ations of lBresolin et al] (|2004l ) and 
IPilvugin fc Thuanl j|2005l ) to determine nebular abundances. 

iModi az et alJ l|2008l ) have investigated the metallicities 
at the sites of type Ic SNe and GRB related SNe and show 
the importance of comparing abundances derived by self- 
consistent methods. By employing consistent abundance in- 
dicators they find that GRB related SNe tend to have signif- 
icantly lower metallicities within their host galaxy environ- 
ments than broad lined type Ic SNe without GRBs. Their 
study illustrates the need to adopt consistent methods and 
compare abundances differentially. 

We chose which stellar tracks to estimate the initial 
masses as follows. We use the o bserved presen t day oxy- 
gen abundances as compiled by iHunter et al.l (|2007h for 
the Sun, LMC and SMC (8.65, 8.35, 8.05) to guide our 
choice of model. For those SNe which have estimated ISM 
oxygen abundances of [O/H] ^ 8.4 we choose to use the 
Z — 0.02 metallicity tracks (solar). For those in the range 
8.2 sC [O/H] < 8.4 we use the Z = 0.008 tracks (LMC). We 
would have used the Z — 0.004 tracks for anything below 
[O/H] < 8.2 (SMC), but none of our targets have such low 
metallicity. 



5 THE MASSES OF THE PROGENITORS OF 
TYPE II-P SUPERNOVAE 

We expect that the progenitors of various subtypes of type II 
SNe are hydrogen rich stars which have evolved from main- 
sequence stars of an approximate initial mass of 8Mq and 
above. If objects are detected in pre-explosion images then 
their colours and luminosities can be determined. If there 
is no star detected at the SN positions then the sensitiv- 
ity of the images can be used to determine an upper lumi- 
nosity limit and hence upper mass limit (see for example 
Smartt et afll2003l : IVan Dvk et~aill2003d ; iMaund fc Smarttl 
20051 ). In the papers presenting the original results, slightly 
different methods of determining the luminosity and mass 
limits from the prediscovery images have been adopted. A 
mixture of 3a and 5a limits have been quoted, uncertain- 
ties treated in varying manners, different stellar evolutionary 
models adopted and distances used which were not always 
the most recent and most accurate. In order to compare the 
sample as a whole, this calls for some homogenisation, and 
we adopt the following method. 

In the cases where there is no detection of a progen- 
itor (13 in total) we determine the upper luminosity limit 
corresponding to the 84 per cent confidence limit. First we 
take the 3a detection limit for each pre-explosion image, 
where this is the detection magnitude in the filter system 



employed. To convert this to a bolometric luminosity one 
requires a measurement of extinction, distance and bolo- 
metric correction (with respect to the filter employed) for 
the progenitor. The la uncertainties of these quantities are 
combined in quadrature to give a total la uncertainty on the 
upper limit. If one assumes that the progenitor star was a red 
supergiant just prior to explosion then the bolometric and 
colou r corrections for an M0 supergiant (jDrilling fc Landold 
l200d ) are appropriate. This assumption seems well justified 
as nearly all the SNe in our sample have been shown to type 
II-P, which require e xtended atmospheres physical ly similar 
to red supergiants (|Arnetd Il980l : IChevalie rf ll976l ). Recent 
detections of the UV shock breakout from two type II-P 
SNe determined the radii of the progenitor stars which adds 
further weight to the idea of the progenitors being red super- 
giants l|Schawinski et al.ll2008l ; iGezari et al.ll2008l ). The un- 
certainty in the bolometric correction is taken to be ±0.3 m 
corresponding to the la range of values for red supergiants 
betw een late-K and late M-type supergiants l|Levesaue et al.l 
2005). Assuming that the uncertainties are representative 
of a normal distribution of measurements, the 84 per cent 
confidence limit for the upper luminosity limit is la above 
the best estimate i.e. there is an 84 per cent chance that the 
progenitor stars have luminosities below this value given the 
individual uncertainties in the calculation. The la distance 
and extinction errors are taken from the quoted sources as 
listed in the notes on the individual events below. These 
84 per cent upper luminosity limits are then plotted on the 
final mass-luminosity plot discussed in Section[3] We deter- 
mine the upper mass limit to be the maximum mass of a 
star which does not have part of its post-He burning track 
within the 84 per cent luminosity limit. 

This method is equivalen t to that previously employed 
by ISmartt et all (|2003l . l2002h and IMaund fc Smarttl (|2005l ) 
(for example) which an exclusion region of the HR diagram 
was determined as a function of effective temperature and 
an upper mass limit from the red supergiant region was de- 
termined. If an /-band (or /-band like) filter is employed 
both methods have the advantage of being fairly insensitive 
to the effective temperature of the assumed red supergiant 
progenitor as the peak of the stellar SED at this temperatur e 
range is -8300A(e.g see Figs. 5 and 6 of lSmartt et~ai1l2003l ). 
In all cases we have revised the distances to the galaxies to 
the most reliable, in our opinion, and most recent in the 
literature. Where no other distance is available we have cal- 
culated a kinematic distance estimate using the host galaxy 
radial velocity corrected for the local group infall into Virgo 
(V v ir from the LEDA0 galaxy catalogue) and a value of the 
Hubble constant of Ho = 72km s _1 Mpc -1 . In such cases we 
employ an uncertainty of the local cosmic thermal velocity 
187kms -1 i|Tonrv et al.ll2000l ). equivalent to ±2.6 Mpc. If 
the value of Ho adopted were either 65 or 85kms -1 Mpc - , 
the systematic effect on the distance scale would provide sys- 
tematic luminosity differences of +0.22 dex and — 0.36dex 
for the five SNe which have kinematic host galaxy distances 
(1999br, 1999ev, 2004dg, 2006bc, 2007aa). We will discuss 
the effects of this systematic difference in Sect .[7] When de- 
termining the extinction in ea ch wavelength band, we use 
the law of lCardelli et alJ (|l989l ). 
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Table 2. The results of the homogeneous reanalysis of the all the SN progenitors. The galaxies, their class, distance and extinction along 
the line of sight to the SNe are listed. The methods employed in the literature to determine distances are listed (TF = Tully Fisher; 
Kin. = Kinematic; Cep. = Cepheid; PNLF = Planetary Nebulae Luminosity Function; TRGB = Tip of the Red Giant Branch; Mean 
= mean of several methods which are detailed in Section 5). Kinematic distances are based on Ho = 72kms _1 Mpc - 1 . The de-projected 
galactocentric radii are calculated as well as the radius with respect to the T2s value (the radius at which the surface brightness drops 
to 25 mag arcsec - 1 ). Oxygen abundances of the galactic ISM at the positions of the SNe are quoted ([O/H] = 12 + log Nq/Nh). The 
final estimated luminosities, or luminosity limits are in solar luminosity units. The zero-age main-sequence (ZAMS) masses and upper 
mass limits as discussed in Sect[5]are listed in the final column. 



Supernova 


SN 


Galaxy 


Galaxy 


Distance 


Ay 


r G 


»"gA25 


[O/H] 


log L/Lq 


ZAMS 




Type 






Class 


Mpc 


ivie Liiou 




(kpc) 




(dex) 


(dex) 


(M Q ) 


1999an 


II 


IC 755 


SBb 


18.5 ± 1.5 


TF 


0.40 ± 0.19 


4.7 


0.82 


8.3 


< 5.16 


< 18 


1999br 


II-P 


NGC 


4900 


SBc 


14.1± 2.6 


Kin. 


0.06 ± 0.06 


3.1 


0.69 


8.4 


< 4.76 


< 15 


1999em 


II-P 


NGC 


1637 


SBc 


11.7 ± 1.0 


Cep. 


0.31 ± 0.16 


1.6 


0.28 


8.6 


< 4.69 


< 15 


1999ev 


II-P 


NGC 


4274 


SBab 


15.1± 2.6 


Kin. 


0.47 ± 0.16 


5.3 


0.46 


8.5 


5.1 ±0.2 


161! 


1999gi 


II-P 


NGC 


3184 


SABc 


10.0±0.8 


Mean 


0.65 ± 0.16 


3.1 


0.30 


8.6 


< 4.64 


< 14 


2001du 


II-P 


NGC 


1365 


SBb 


18.3 ± 1.2 


Cep. 


0.53 ± 0.28 


14.7 


0.53 


8.5 


< 4.71 


< 15 


2002hh 


II-P 


NGC 


6946 


SABc 


5.9 ±0.4 


Mean 


5.2 ± 0.2 


4.1 


0.45 


8.5 


< 5.10 


< 18 


2003gd 


II-P 


NGC 


628 


Sc 


9.3 ± 1.8 


Mean 


0.43 ± 0.19 


7.5 


0.58 


8.4 


4.3 ±0.3 


7+6 
'-2 


2003ie 


II? 


NGC 


4051 


SABb 


15.5 ± 1.2 


TF 


0.04 


7.3 


0.66 


8.4 


< 5.40 


< 25 


2004A 


II-P 


NGC 


6207 


Sc 


20.3 ± 3.4 


Mean 


0.19 ± 0.09 


6.7 


0.79 


8.3 


4.5 ±0.25 


7+6 
'-2 


2004am 


II-P 


NGC 


3034 


Sd 


3.3 ±0.3 


Cep. 


3.7 ± 2.0 


0.64 


0.14 


8.7 


Cluster 


X2±l 


2004dg 


II-P 


NGC 


5806 


SBb 


20.0 ±2.6 


Kin. 


0.74 ± 0.09 


4.3 


0.50 


8.5 


< 4.45 


< 12 


2004dj 


II-P 


NGC 


2403 


SABc 


3.3 ±0.3 


Cep. 


0.53 ± 0.06 


3.5 


0.37 


8.4 


Cluster 


15 ±3 


2004et 


II-P 


NGC 


6946 


SABc 


5.9 ±0.4 


Mean 


1.3 ± 0.2 


8.4 


0.92 


8.3 


4.6 ±0.1 


9±f 


2005cs 


II-P 


NGC 


5194 


Sbc 


8.4 ± 1.0 


PNLF 


0.43 ± 0.06 


2.7 


0.22 


8.7 


4.25 ±0.25 


7+3 
' -1 


2006bc 


II-P 


NGC 


2397 


SBb 


14.7 ±2.6 


Kin. 


0.64 


1.4 


0.30 


8.5 


< 4.43 


< 12 


2006my 


II-P 


NGC 


4651 


Sc 


22.3 ±2.6 


TF 


0.08 


4.4 


0.37 


8.7 


< 4.51 


< 13 


2006ov 


II-P 


NGC 


4303 


SBbc 


12.6 ±2.4 


TF 


0.07 


2.3 


0.26 


8.9 


< 4.29 


< 10 


2007aa 


II-P 


NGC 


4030 


Sbc 


20.5 ±2.6 


Kin. 


0.09 


10.3 


0.91 


8.4 


< 4.53 


< 12 


2008bk 


II-P 


NGC 


7793 


Scd 


3.9 ±0.5 


TRGB 


1.0 ± 0.5 


3.9 


0.66 


8.4 


4.6 ±0.1 


9±t 



In the cases where we have a direct detection of the 
progenitor (4 in total) the uncertainties in the luminosity 
are trivially determined and discussions of the individual 
events are listed below. 



Two others (2004dj and 2004am) fall on bright, com- 
pact star clusters which are not resolved into individual stars 
llMafz-Apellaniz et al.ll2004l ; IWang et al1l2005l : IVinko et all 
|200S| , Mattila et al. 2008, in prep). These papers have de- 
termined the total cluster mass and age and hence the turn 
off mass at the top of the main- sequence. From this the 
mass of the progenitor has been determined. Hence the stel- 
lar mass determination is somewhat indirect and relies on 
the assumption that the timescale of star formation in the 
cluster is significantly less than the current estimated age. 
The results are based on population synthesis codes which 
use different individual stellar evolutionary codes as input t o 
those we have emp loyed. The iMaiz-Apellaniz et "all (|2004l ); 
IVinko et al.l (|2008T ) results are based on the stellar synthe- 
sis codes starburst 99 which uses the Geneva models (e.g. 
ISchaller et al.| [l992) as input, and the discussion in Sec- 
tion l6.ll indicates that the choice of stellar model does not 
introduce significant uncertainties. Hence although the anal- 
ysis method, and hence mass determination, is different for 
these two events, we believe the mass estimates are worth 
including in this compilation. If they are left out, the main 
conclusions of this paper are not altered in any significant 
way. 



5.1 1999an 

IC755 is an SBb spiral with M B = -18.85 (from LEDA). 
As no direct abundance study of this galaxy has been done, 
we attempt to infer a probable abundance at the position of 
the pr ogenitor from the O/H— Mb relation of lPilvugin et al.l 
(2004). For galaxies between —20 < M B < -19, the charac- 
teristic oxygen abundance (the oxygen abundance at a galac- 
tocentric distance of r = 0.4r2s) is typically in the range 
8.5 ± 0.2 dex. The mean abundance gradient of this sample 
is —0.5 ± 0.3dex/r25. Hence at a de-projected distance of 
0.8r25, the metallicity of the progenitor star of SN1999an 
can be approximated at 8.3 dex. This is somewhat uncer- 
tain given the large uncertainties on the gradient and the 
range of characteristic abundances and the error is likely to 
be ±0.3 dex. However it is the best estimate that can be 
derived with the current data in the literature. 

There is no detection of a prog enitor star in the WFPC2 
pre- explosion images prese nted by iMaund fc Smarttl (|2005l ) 
and IVan Dvk et"atl l|2003al ). Both studies calculated similar 
se nsitivity limits for the im ages, and we adopt the 3<r limit 
of lMaund fc Smarttl (2005) of m,F606w = 24.7. Solanes et al. 
(2002) report a mean distance modulus for the host galaxy 
IC755 of fi = 31.33 ± 0.18 or d = 18.45 ± 1.5Mpc. Apply- 
ing a line of sight extin ction of E(B - V) = 0.13 ± 0.06 
(jMaund fc Smarttl 120051 ) and assuming an M-type super- 
giant as the progenitor (for the colour correction between 
F606W and Johnson V ; see Maund & Smartt 2005) results 
in absolute upper limit of My = —6.46 ± 0.26. For an M0 
supergiant this corresponds to an upper luminosity limit of 
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logL/L0= 5.00 ± 0.16, and an 84 per cent confidence limit 
of logL/L = 5.16. From Fig. [T]this implies an upper mass 
limit of 18M . 



5.2 1999br 

NGC4900 is an SBc spiral, and similar to the case of IC755 
discussed above it does not have a published abundance 
study. The same arguments as in Section l5TTl can be used 
(NGC4900 has Mb = -19.05) to infer an oxygen abun- 
dance at the galactocentric radius of SN1999br (0.69r2s) of 
approximately 8.4 dex. 

Wit h no detection of a pro genitor object at the SN 
position, iMaund fc Smarttl (|2005l ) place an upper limit on 
the magnitude of any progenitor o f mF606w = 24 . 9. This 
is significantly brighter than that of I Van Dvk et all (|2003al ) 
wiF606W = 25.4, and we adopt the former as the more conser- 
vative result. For the host galaxy (NGC4900) only a kine- 
matic distance modulus is available, with the Virgo infall 
corrected velocity giving d = 14.1 ± 2.6 Mpc. The extinc- 
tion to this event appears very low (IMaund fc Smarttll2005l; 
I Van Dvk et al.ll2003al ; iPastorello et al.ll2004l ) and we adopt 
the foreground value quoted in these papers of E(B — V) = 
0.02 ±0.02. As for SN1999an, we assume that the progenitor 
was a red supergiant and apply a colour correction and bolo- 
metric correction to determine an upper luminosity limit of 
logL/L0= 4.55 ± 0.21 and an 84 per cent confidence limit 
of logL/L0= 4.76. From Fig[T]this implies an upper mass 
limit of 15M Q . 



5.3 1999em 



Ivan Zee et alj l| 19981 ) have published line strength measure- 
ments and abundance s for 15 Hu reg i ons in NGC1637. Us- 
ing the calibration of iBresolin et al.1 l|2004t ) we have rede- 
termined the abundance gradient and at the galactocentric 
distance of SN1999em, the metallicity is 8.6 ± 0.1 dex. The 
nearest Hu regions to 1999em are 510 pc and 794 pc, when 
de-projected, from the site of SN1999em and have oxygen 
abundances of 8.5 and 8.7 dex respectively. Hence we adopt 
8.6 dex. 

The updated distance to NGC1637 of d = 11.7 ± 
1 Mpc is taken from the Cepheid variable star estimate of 
iLeonard et~alj (|2003h and the reddening value of E(B-V) = 
0.1 ± .05 is adopted from lBaron et alTl|2000l ). ISmartt et al.l 
(2002) present deep ground based images of NGC1637 be- 
fore explosion in VRI filters and from these results we have 
determined a 3a upper limit of I — 23. The /-band is the 
most sensitive to red supergiant progenitors and between 
the supergiant spectral types of K2-M4 this corresponds to 
an upper luminosity limit of logL/Lo= 4.54 ± 0.15, and an 
84 per cent confidence limit of log L/Lq— 4.69. From Fig[T] 
this implies an upper mass limit of 15M0. 



5.4 1999ev 

NGC4274 is an SBab spiral, and also has no abundance 
study of its Hu regions. The same arguments as in Sec- 
tional] can be used (NGC4274 has M B = -20.18) to in- 
fer an oxygen abundance at the galactocentric radius of 
SN1999ev (0.46r 25 ) of approximately 8.5 dex. 



SN1999ev was recovere d in late, deep HST ACS images 
bv lMaund fc Smarttl <|2005l ) and is coincident with a progen- 
itor object fou nd on a pre-explosion WFPC2 F555W im- 
age. Although IVan Dvk et alj (l2003al ) originally suggested 
two other stars as possible progenitors, the HST follow- 
up clearly ruled this out and points to the object of 
magnitude mFsssw = 24.64 ± 0.17, at 4.8cr significance 
l|Maund fc Sma rtt 2005). There is no distance measurement 
to the galaxy NGC4274 apart from a kinematic estimate, 
which is d = 15.14 ± 2.6 Mpc, fro m LEDA (Virgo infall cor- 
rected) . [Maund~&^Smartl] |2005l ) determined the extinction 
to the nearby stellar population of E(B — V) = 0.15 ± 0.05. 
We again assume that the progenitor was a red supergiant 
and apply a BC of — 1.3±0.3 to determine a final luminosity 
of logL/L0= 5.1 ± 0.2. The tracks in Fig[T] imply the star 
would have been of mass 1615m 



0- 



5.5 1999gi 

ISmartt et al.l J200ll) sugg ested that the Hu region number 
3 of lZaritskv et al.l (|l994 ) at a position of 68"N 0"E is coin- 
cident with the star-forming region, o r OB association that 
hosted SN1999gi. The calibration of I Bresolin et al.1 (j2004l ) 
using the R23 value of IZaritskv et aL ( 19941 ) gives an oxy- 
gen abundance of 8.6 dex. 

A study of the progen itor site of SN1999gi was carried 
out bv lSmartt et al.l (|200ll ). but the distance to this galaxy 
was th en im proved i n a com pilation study of lLeonard et alj 
(|2002|) and IHendrvl (|2006r ). Here we adopt the result in 
Hcndrv (2006) which is a mean of four estimates d = 10.0 ± 
0.8 Mpc, and the extinction of lLeonard et al] E(B — V) = 
0.21 ± 0.05. The four distance methods detailed in IHendrvl 
(2006) are Tully Fisher, expanding photosphere method 
(EPM), kinematic an d the tertiary distance indicators of 
Ide Vaucouleursl d 19791). The 3<r detection limit determined 
by ISmartt et a l. ( 2001) is 771F606W = 24.9, which results 
in an upper luminosity limit for an M-type supergiant of 
logL/L©^ 4.49 ± 0.15, after the colour and bolometric cor- 
rections are applied. This gives an 84 per cent confidence 
limit of logL/L0= 4.64 and, from Fig. [TJ this implies an 
upper mass limit of 14M© . 



5.6 2001du 



As discussed in ISmartt et al.l (|2003l ) the Hu region RW21 
l|Rov fc Walshlll997l ) is 1.5" from SN2001du, and it is likely 
the metallicity of this region is representati ve of the pro- 
enito r star composition. T he calibration of Bresolin et al.l 
20041 ) to the R23 value of iRov fc Walsh! j|l997l ) gives an 
oxygen abundance of 8.5 dex. 

The host galaxy NGC1367 was observed as part of 
the HST Cepheid Key Project, hence the most accurate 
and r ecent distance estimate is taken from IPaturel et al.1 
(2002), fi = 31.31 ± 0.15. The extinction to ward the SN 
was m easured by three different methods by ISmartt et al.1 
(|2003h to be E{B - V) = 0.17 ± 0.09, giving A 7 = 
0.25 ± 0.13, which is similar to that adop t ed (A i ~ 
0.2) by IVan Dvk et all (|2003d ). ISmartt et all (|2003l ) and 
IVan Dvk et al. (2003c) presented pre-explosion images in 
the WFPC2 filters F336W, F555W, F814W and the most 
sensitive of these to red supergiants is the F814W. We de- 
termine the 3ct upper limit from the Smartt et al. results 
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to be I = 24, similar to the sensitivity mFsi4,w = 24.25 
of I Van Dvk et al.| (|2003d V Between the supergiant spectral 
types of K2-M4 this corresponds to an upper luminosity 
limit of logL/L0= 4.57 ± 0.14, and an 84 per cent confi- 
dence limit of logi/L0= 4.71. From Fig[T]this implies an 
upper mass limit of 15M©. 



5.7 2002hh 

None of the 9 H i i regions in NGC6946 compiled by 
iPilvugin et"al]||2004l ) are near the location of 2002hh. Hence 
we use the abundance gradient determined bv lPilvugin et all 
( 2004) and the de-projected galactocentric radius of the SN 
position to determine t he likely metal l icity. As discussed 
above, the calibration of lPilvugin et al.l 1120041 ) is sim ilar to 
the simple linear calibration of lBresolin et alJ (2004), hence 
these should be on a similar scale. We determine an oxygen 
abundance of 8.5 dex. 

A deep pre-explosion i'-band archive image of NGC6946 
from the Isaac Newton Telescope Wide Field Camera (INT- 
WFC) will be presented in a forthcoming paper (See 
Sect l5.14]) . Although this SN suffered significant extinction, 
the proximity of the galaxy and the depth of the 3600s i'- 
band image still places useful restrictions on the progenitor 
star. The 3600s image is composed of 6 x 600s exposures, 
with a final image quality of 1". There is no object visible 
at the position of SN2002hh, and the 5<r detection limit for 
a point source was estimated to be iccd=22.8. This instru- 
mental magnitude can be converted to a standard I using 
the well calibrated co lour transformations for the INT-WFC 
ijlrwin fc Lewis! 1200 1 ) Q We e mploy the reddening law de- 
termined in iPozzo et al.l J2006) to estimate the extinction in 
the J-bandof Ar = 2. 1±0. 3. A distance of d = 5.9±0.4Mpc 
is used which is a mean of the distance values from the com- 
pilation of Botticella et al. (2009 in prep.) using the methods 
of Tully Fisher, brightest supergiants, sosies, PLNF, EMP 
(applied to 1980K) and standard candle method (SCM) ap- 
plied to SN2004et. SN2002hh appears to be a no rmal II- P, 
but behind a large dust pocket (|Pozzo et al . 2006), hence we 
assume the progenitor was a red supergiant of type between 
K0-M5. The falling bolometric correction combined with the 
rising intrinsic (V — I) between K0-M5 means that the bolo- 
metric luminosity limit stays approximately constant in this 
spectral range at log L/Lq = 4.9 ± 0.2. Hence the 84 per 
cent confidence limit is logL/L0= 5.1 and from Fig[T]this 
implies an upper mass limit of I8M0. We note that this is 
co nsistent with the p rogenitor mass of I6-I8M0 estimated 
bv lPozzo et "all (|2006T ) from the [0 1] AA6300,6364A doublet. 



5.8 2003gd 

None of the previously catalogued H n regions in NGC628 
which have spectra and the R23 ratio measured are par- 
ticularly near the spatial position of SN200 3gd. Hence we 
use th e abundance gradient determined by IPilvugin et al.l 
|2004h and the de-projected galactocentric radius of the SN 
position to determine the metallicity at this position. The 
parameters derived are Tabled with an abundance of 8.4 dex 
derived. 



http: / /www. ast.cam.ac.uk/ ~wfcsur/index.php 



The p rogenitor star was de tected by ISmartt et al.l 

l|2004h and IVan Dvk et all (|2003d ). and an extensive com- 
pilation of distance measurements to NGC628 and r edden - 
ing towards the SN was carried out in iHendrv et all (|2005l ). 
Those distance and e xtinction value s deter mined were close 
to those employed in ISmartt et al.l (|2004|) to estim a te the 
progenitor luminosity and mass. The IHendrv et all (|2005h 
distance listed in Table[2] is a mean of the three meth- 
ods : kinematic, brightest supergiants and standard candle 
method (applied to SN2003gd). The intrinsic (V - I) = 
2.3±0.3 colour is consisten t with a supergiant in the spectral 
type range K5-M3, which ISmartt et al.1 (|2004l ) used to de- 
termine a luminosity of log L/L0= 4.3 ± 0.3. In the diagram 
of Fig[T] the best value of 4.3 is closest to the termination 
point of the 7M0 track, if we assume the progenitor did not 
go through 2nd dredge-up. The uncertainties would bracket 
the 5M0 and 13M0 post-He burning tracks, hence we adopt 
the value 7t|M0. Although the most likely value is below 
the lowest mass that is normally assum ed possible to pro- 
vide an iron core- collap se (8-lOM ; see iHeger et~aTI 120031 ; 
lEldridge fc Toutl [2004b), the range of masses comfortably 
brackets the theoretically predicted limits. 

5.9 2003ie 

NGC4051 is a Seyfert 1 galaxy of morphological type SABb 
and has no published abundance study of its Hn regions. 
Hence we can estimate a probable abundance at the po- 
sition of the progenitor as in Sect . 15.11 The galaxy has a 
Mb = —20.3, using our adopted distance and the corrected 
B— band magnitude given in LEDA. At this magnitude 
the the characteristic oxygen abundance (at a galactocen- 
tric distance of r = 0.4r2s) is approximately 8.5 ± 0.2 dex 
l|Pilvugin et al. I bin. || i Again using the typical gradient of 
the galaxy sample (as in Sect. l5.ll of —0.5 ± 0.3 dex/r25, the 
oxygen abundance at 0.66r25 is approximately 8.4 dex. As 
stated above this is quite uncertain (±0.3 dex) given the 
lack of detailed study of the galaxy but does show that it 
is unlikely to be a particularly low metallicity environment. 
The SN was not studied in great detail by any group (as 
far as we know), but a single photospheric spect rum shows 
P-Cygni features of H 1 (|Harutvunvan et al.ll2008l). The best 
match for the spectrum found bv lHarutvunvan et al.l (|2008l ) 
is that of 1998A, which itself appears like a 1987A-type 
event. Hence this event may not be a normal II-P SN and 
we have no lightcurve information to consider. As we shall 
see below, the mass limits for the progenitor are not partic- 
ularly restrictive and if the object were to be left out the 
conclusions of the paper would be unchanged. 

The pre-explosion site of SN 2003ie was recovered in 
archive r'-band observations of NGC4051 taken with the 
INT-WFC. This image was taken on 1999 November 11, 
with an exposure time of 900s and image quality of 1.1 arc- 
sec. We determined the position of SN 2003ie within an error 
circle of 0.17 arcsec (using an image of the SN from 2003 pro- 
vided to us by Martin Mobberley). There is no progenitor 
object detected within this error circle and we derive a 3<r 
detection limit of r' = 23. This instrumental magnitude was 
converted to a standard R using th e well calibrated colo ur 
transformations for the INT-WFC l| Irwin fc Lewis! l200ll ) 7 . 
giving R = 22.65. Pierce & Tully (1988) calculate the dis- 
tance to the Ursa Major Cluster to be 15.5 ± 1.2 Mpc us- 
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ing the Tully-Fisher method and we adopt this distance for 
NGC4051. We have no measure of the internal extinction 
towards this SN and simp ly adopt the Galact ic extinction 
value of E(B-V) = 0.013 jSchlegel et al.lll998l ). We assume 
once again that the progenitor is a red supergiant and apply 
appropriate bolometric and colour corrections to determine 
a luminosity limit of log L/Lq = 5.26 ± 0.14 and an 84 per 
cent confidence limit of logL/L© = 5.40. From Fig. [1] this 
implies an upper mass limit of 24M©. 

5.10 2004A 

There are no measurements of H II regions in NGC6207 and 
we use the arguments presented in iHendrv et al.l l)2006l l to 
estimate the metallicity at the galactocentric distance of 
SN2004A. This paper based the results on ty pical abundance 
gradients measured bv lPilvugin et all lj2004T ) hence a sain the 
estimate is on the same scale as the rest of the values. We 
point out that there is a typographical error in that paper, 
where 7?25 is quoted as 4kpc, whereas it should be 8.6 kpc. 
However, repeating the same method, this does not change 
the calculated metallicity at the position of SN2004A which 
we estimate as 8.3 dex. The IHendrv et ail (2006) distance 
listed in Tabled is a mean of the three methods : kinematic, 
brightest supergiants and standard candle method (applied 
to SN2004A). 

A faint object is d etected at the position of SN2004A 
in IHendrv et alj (|2006T ). claimed as a 4.7<j detection in the 
F814W filter, and it is not detected in the F435W or F555W. 
If we assume this detection to be valid, it provides a blue 
limit for the colour of the progenitor and hence a star in the 
spectral range G5-M5. This gives a bolometric luminosity 
in the range logL/L/0= 4.5 ± 0.25. On Fig[T]this implies a 
best estimate of 7Mq , and the errors bracket the 5Mq and 
13Mq post-He burning tracks. Hence we adopt 7^}Aq. If 
the detection is not valid then the I-band detection sensitiv- 
ity implies an 84 per cent confidence limit of log L/Lq=4.75 
and an upper limit of 13M Q . 



5.11 2004am 

There is no extensive published study of SN2004am to 
date, which is surprising given its proximity and the fact 
it is the only optically discovered SN in the starburst M82 
(NGC3034). However it is clearly a II-P from the unfiltered 
magnitudes of lSinger et al.l d2004h which stay co nstant for 76 
days, and the spectrum of lMattila et al ] <|2004h . Alignment 
of post-explosion near-IR images and HST pre-explosion 
images shows SN2004am is spatially coincident with the 
well studied super star cluster M82-L (Mattila et al., in 
prep). The distance to M82 is assumed to be that of the 
M81 group, estimate d from Cepheids in M81 (NGC3031; 
iFreedman et al.ll200ll ). 

A new study of M 82-L has recently been carried out 
by lLancon et all (|2008T l who modelled the integrated near 
infra-red 0.8 — 2.4/jm spectra. They used a population syn- 
thesis code (PEGASE.2) and a new library of red super- 
giant observational and theoretical spectra to determine the 
age of M82-L. The fit to the overall SED and the individ- 
ual molecular absorption features is impressive and gives 
an age estimate of 18lg 7 Myrs. With a lower than normal 



value of Ry ~ 2.4 — 2.7, lLancon et all l|2008h can also re- 
produce the optical SED of the cluster down to 6000A (the 
Ay value in Tabled is taken from Lancon et al.). This age 
is somew hat younger than 651 35 Myrs that was first in- 
ferred bv lSmith et ail |2006) using only a limited range op- 
tical spectrum and th e spectral synthesis c ode starburst99 
l|Leitherer et all 1 19991) . lLancon et ail (|2008h point out that 
by using a low value of Rv they can reconcile the optical 
SED and the NIR molecular bands with the younger age 
and their updated spectral modelling technique and firmly 
exclude an age of 60Myr. The cluster age provides quite a 
strong constraint on the mass of the progenitor star, assum- 
ing that the cluster formed coevally and the progenitor's 
age is similar to that of the cluster. The STARS models pre- 
dict that the cluster ages correspond to lifetimes of stars of 
masses 12 l 7 Mp. The models u sed in the pop ulation synthe- 
sis cod e of lLancon et al.1 l|2008l ) were those of iBressan et al.l 
i| 19931 ). which give very similar age-mass relationships to the 
STARS code (the uncertainty on the derived mass due to 
choice of code is within the error range). 

In all of the above we have assumed solar metallicity for 
the stellar evolution ary tracks is appropriate. The fitting of 
lLancon et al.l (|2008l ) implies that this is appropriate. Also 
two recent papers have speculated on the abundances in the 
nuclear regions of M82, in environments close to su per star 
cluster M82-L (|Smith et al.ll2006l ; lQriglia et alj|2004 ). There 
is some uncertainty and difference in the abundances derived 
but the stell ar abundances of red supergiants in these inner 
regions from lOriglia et al.l l|2004l) suggest a solar like oxy- 
gen abundance. The photospheric abundance in such objects 
are likely to be applicable to M82-L and the progenitor of 
SN2004am. Hence in the age estimations we chose the tracks 
close to 8.7 dex as the most appropriate. 

5.12 2004dg 

NGC5806 is an SBb spiral, and also has no abundance study 
of its Hn regions. The same arguments as in Section [5Tl can 
be used (NGC5806 has Mb = -19.86) to infer an oxy- 
gen abundance at the galactocentric radius of SN2004dg 
(0.56r25) of approximately 8.5 dex. 

The pre-explosion site of SN 2004dg was imaged using 
both WFPC2 (2001 July 5) and ACS (2004 April 3) cameras 
on-board HST (the SN was discovered on 2004 July 31). The 
WFPC2 exposure times were 460s in F450W and F814W fil- 
ters. For the ACS images had total exposure times of 700s 
in F658N and 120s in F814W. We re-observed SN2004dg on 
(2005 March 10) with the ACS camera (in F435W, F555W 
and F814W, as part of GO10187) and recovered the SN at 
transformed magnitudes of B = 22.2, V = 20.8, / = 19.1. 
Alignment of the two sets of images allowed us to locate the 
position of the S N on the pre-explo s ion images to within 
0.015 arcsec (see iMaund et ail 120051 ; ICrockett et al.1 120081 . 
for details of alignment procedures) Within this error cir- 
cle there was no detection of a progenitor star in any of the 
image and filter combinations. A progenitor star was not de- 
tected at the SN position, therefore a 3a detection limit of 
mF8i4,w = 25.0 was determined. There is no distance mea- 
surement to the galaxy apart from a kinematic estimate, 
which is d = 20.0 ± 2.6 Mpc, from LED A (Virgo infall cor- 
rected). The total reddening towards SN 2004dg was esti- 
mated to beE(B-V) = 0.24±0.03, giving Ai = 0.36±0.04. 
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Again we assume a red supergiant progenitor and find a lu- 
minosity limit of log L/Lq = 4.28 ±0.17 and an 84 per cent 
confidence limit of log L/Lq = 4.45. From Fig.[T]this implies 
an upper mass limit of 12Mq. 



5.13 2004dj 

As discussed in Sect.0 SN2004dj fell on the com pact 
star cluster i denti fied by IMaiz-Apellaniz et alj l|2004h and 
IWang et all ([2005) and the analysis used to determine a 
progenitor mass is different to the direct identification and 
direc t upper luminosity limits for the other SNe presented 
here. IMaiz-Apellaniz et al.l (|2004| ) determine an age of the 
compact star cluster of 14Myrs a nd hence a main-se quence 
mass of 15M0 for the progenitor. IWang et all (|2005l ) deter- 
mine an age of around 20Myrs and hence a m ain-sequence 
mass of 12M©. A new and improved study by IVinko et al.l 
(2008) using new UV observations of the cluster and ex- 
tensive comparison of SEDs based on different model at- 
mospheres and evolutionary tracks suggests a most likely 
turn off mass (and hence progenitor ZAMS mass of between 
12-2OM . This age (10-16 Myrs) is consistent with the lack 
of Ha emis sion seen in the cluster spectrum of IVinko et al.l 
l|200fil) and [Humphreys fc Aaronsonl (119871 ). as the ionizing 
O-stars have died out. Hence we favour the older age and 
will adopt 15 ± 3M0 as the progenitor ma ss (if we adopt 
15^3^40 as suggested by I Vinko et al.l [ 20081 , it does not af- 
fect any of the results below). The distance to NGC2403 
of 3.3 ± 0.3 Mpc is f rom the HST Cepheid Key Project 
llFreedman et al.ll200ll) 



Mai'z-Apellaniz et alj (|2004l ) adopted solar abundances 



in using the Geneva tracks of starburst99. On closer inspec- 
tion this may b e too high. We used t he abundance gradient 
determined by iPilvugin et al.l (|2004l ) and the de-projected 
galactocentric radius of the SN position to determine the 
metallicity at this position of 8.4 dex. Although the Hn re- 
gions in this galaxy have been studied extensively the only 
region which is physical ly close to the posit ion of 2004dj is 
that of VS44 studied bv lGarnett et al.l (|l997f) . and even that 
is around 600pc from t he cluster that hosted SN2004dj. The 
-R23 ratio provided by iGarnett et al.l (|l997l ) also gives an 
abund ance of 8.4 dex with the calibration of iBresolin et al.1 
(2004), in good agreement with the abundance gradient 
measuremen t. Although our oxygen abun dance is below that 
employed bv lMafz-Apellaniz et al.1 l|2004h . this does not sig- 
nificantly affect the age (and hence turn-off mass) estimate 
when we compare STARS models of s uch different metallic- 
ities. We note that IVinko et al.l (|2008h favour a solar metal- 
licity in their SED fits. 



5.14 2004et 

As discussed for SN2002hh C Sect. RTF) . there is no Hn re- 
gion near the galactic position of 2004et, which is some way 
from the centre of NGC6946. We use the same method as 
for 2002hh to determine a metallicity typical for the galac- 
tocentric radius of 2004et of 8.3 dex. The adopted distance 
to NGC6946 of d = 5.9 ± 0.4 Mpc is discussed in SectlBTTl 

iLi et al.l l|2005h presented the detection of a candi- 
date progenitor star of SN2004et in ground-based images 
from the Canada France Hawaii Telescope (CFHT) in both 



UBVR and u'g'r' filters. They suggested it was a yellow 
supergiant as the BVR colours were matched with a G- 
ty pe supergiant SEP. From t hese colours an d the models 
of lLeieune fc Schaererl l|200ll) . ILi et al.1 l|2005h determine a 
mass in the range I5I2M0 . 

How ever it is now clear that the putative source de- 
tected bv lLi et"ail (|2005l) was not a sing l e star . Adaptive op- 
tics images of the site bv lCrockett et all (|2009f l using Gemini 
North show that the source breaks up into several stars. In 
addition BVR images taken 3 years after explosion show 
identical colours to the pre-explosion object, indicating that 
the progenitor star was not detected in the pre-explosion 
frame. A deep i'-band image (the same image as discussed 
above for SN2002hh) does show a clear detect ion of a progen- 
itor co mpared to the late time /-band image. ICrockett et al.l 
(2009) use the i'-band detection (after converting to John- 
son; I — 22.06 ±0.12) and limits on the BVR magnitudes to 
infer that the progenitor was a red supergiant with (R — I) 
> 1.80±0.22. This implies an M4 spectral type or later, giv- 
ing a bolometric magnitude of of Afboi = —6.73 ±0.22 and a 
progenitor luminosity of logL/L = 4.59±0.09. Comparing 
this to stellar models of LMC metallicity, we estima te its 
initial mass to be 9~t\ Mq (as in lCrockett et af1l2009l ). 



5.15 2005cs 



iMaund et all l|2005h have estimated the abundance at the 
galactocentric radius of SN2005cs in an identical man- 
ner as we have employed consistentl y in this pape r , usin g 
the NGC5194 abundance gradient of lBresolin et al.l l|2004h . 
hence it is already on our common calibration scale. They 
determine 8.66 ± 0.11 dex, which we adopt in this paper. 

The detect ion of the pro g enitor of SN2005cs is well 
documented by lEldridge et all <|2007h . ILi et all (|2006h and 
IMaund et al. I il2005T l which all give s imilar mass estimates in 
the range 7-lOM . lEldridge et al.1 (120071 ) have recently re- 
analysed all of the available photometry from the initial two 
discovery papers and suggested that the progenitor could 
not have been a super-AGB star that has gone through 
2nd dredge up. This analysis was done with the STARS 
code in an identical manner as this study and they found 
likely progenitor r ange of 6 — 8M0 (ass uming a distance of 
8.4±1.0Mpc). The lMaund et all l|2005l ) luminosity estimate 
for the progenitor is log L/Lq = 4.25 ± 0.25 and using the 
STARS tracks employed here this would suggest a mass of 
7l 1M0 . If the closer distance of 7.1±1.2 Mpc is chosen (from 
the mean of the compilation of iTakats fc Vinkdl2006h then 
the best estimate of mass would reduce slightly to around 
6M©. This would be rather low, but the uncertainty on the 
upper bound (3M©) would still place it comfortably within 
the normal theoretical ranges for core-collapse. 



5.16 2006bc 

NGC2397 is an SBb spiral, and is the final galaxy in this 
sample which has no published abundance study of its Hn 
regions. The same arguments as in Section [5~T1 are employed 
(NGC2397 has Mb = —19.67) to infer an oxygen abundance 
at the galactocentric radius of SN2006bc (0.3r2s) of approx- 
imately 8.5 dex. 

The pre-explosion site of SN 2006bc was imaged using 
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WFPC2 (2001 November 17) on-board HST with exposure 
times of 460s in each of the F450W and F814W filters (the 
SN position fortunately fell on the PCI chip). We re ob- 
served SN2006bc on (2006 October 14, as part of GO 10498) 
with the ACS Wide Field Camera (WFC) in three filters 
F435W 1400s, F555W 1500s and F814W 1600s). Aligning 
the before and after explosion images allowed us to locate 
the position of the SN o n the pre-explosion images to within 
0.024 arcsec (again see iMaund et al.l 120051 : ICrockett et al.l 
120081 . for details of alignment procedures). Within this error 
circle there was no detection of a progenitor star in any of 
the image and filter combinations. At the progenitor posi- 
tion we determined a 3er detection limit of mFsuw = 24.45. 
There is no distance measurement to NGC2397 apart from 
a kinematic estimate, which is d — 14.7 ± 2.6 Mpc, from 
LEDA (Virgo infall corrected). The Galactic extinction i s 
estimated to be E(B - V) = 0.205 (|Schlegel et all 1 19981 ). 
Assuming that the progenitor was a red supergiant we find 
a luminosity limit of log L/Lq = 4.23 ± 0.20 and an 84 per 
cent confidence limit of logL/L© = 4.43. From Fig. [1] this 
implies an upper mass limit of 12Mq. 

5.17 2006my 

SN2006my occurred at tg/t^s = 0.37 almost exactly at the 
galactocentric radius of th e characteristic oxyge n abundance 
measured in NGC4651 bv lPilvugin et al.l l|2004h . They mea- 
sure a_rahMof_8 : 7dex at this position. 

iLi et al.l l|2007h claim the detection of a red supergiant 
progenitor of SN 2006my in pre-explosion HST/WFPC2 ob- 
servations of NGC4651. In order to determ ine the position 
of the SN on the pre-explosion HST images ILi et all (|2007l ) 
aligned these images with ground based observations of 
the SN from the Canada- France-Hawaii Telescope (CFHT). 
They derived an initial mass of M = 10±^M© for the object 
which they find coincident with the S N position. 

However in an improved analysis, iLeonard et al.l (|2008T ) 
have shown that this is unlikely to be correct and the pro- 
genitor star is most likely not detected in the pre-explosion 
images. They used HST images of much higher resolution 
than CFHT which allows for object positions to be more 
accurately measured and ultimately leads to a more reli- 
able transformation between the coordinate systems of the 
pre- and post-explosion images. They find that the offset be- 
tween the SN and possible progenitor position is too large 
to support the claim that the two objects are associated (at 
about the 96% confidence level). In a completely indepen- 
dent manner, we used similar data to carry out the same im- 
age alignment and the det ails of this analysis a re presented 
in our companion paper (|Crockett et al.l l2009h . Using the 
HST post-explosion to HST pre-explosion transfo rmation, 
we als o find that the progenitor object proposed bv lLi et"al] 
(|2007h is ~74 mas from the transformed SN position. Given 
our total astrometric error this is approximately a 1.8<r sep- 
aration. Hence we also find that this object is unlikely to be 
the progenitor of SN 2006my. Most likely it is not and the 
progenitor is undetected in the images, so we derive a 3<r 
detection limit of mF8i4w = 24.8. 

8 This constituted an ESA Photo Release : 
http://www.spacetelescope.org/news/html/heic0808.html 



Solanes et al. (2002) have collected Tully-Fisher dis- 
tance estimates for NGC4651 from seven different sources 
and derive a mean distanc e modulus fi = 31.74 ± 0.25 (or 
d = 22.3 ± 2.6 Mpc). As in lLi et all ^OPTf ) we apply only a 
correction for the Galactic extinction of E(B — V) = 0.027. 
Assuming that the progenitor star was a red supergiant we 
derive a luminosity limit of log L/Lq = 4.35 ± 0.16 and an 
84 per cent confidence limit of log L/Lq — 4.51. From Fig.[T] 
this implies an upper mass limit of 13M Q . 



5.18 2006ov 

iPilvugin et ail (|2004T ) redetermined the oxygen abundance 
gradient, and our calibrations are on this equivalent scale, 
hence we use the galactocentric radius 0.26r25 and abun- 
dance gradient of Pilyugin et al. to determine an oxygen 
abundan ce o f 8.9 d ex. which is the highest in this sample. 

ILi et alj l|2007h report the detection of a red supergiant 
progenitor of Mzams = I5I3M0 for SN 2006ov in archival 
HST/WFPC2 observations of NGC4303. In this case the 
pre-explosion frames were aligned with HST observations of 
the SN in order to pinpoint the position of the progenitor 
on the archival images. Having perform ed PSF-fitting pho- 
tometry using HSTphot (|Dolphuj|2000h without detecting a 
progenitor star, it was noticed by Li et al. that a significant 
point source was stil l visibl e in the residual image close to 
the SN site. ILi et alj l|2007l ) suggest that this object in the 
residual image is in fact coincident with the SN position, and 
claim that by forcing HSTphot to fit a PSF at this position 
they detect an object of 6.1<r significance in the F814W and 
F450W observations. 

We have repeated the alignment of the pre- and post- 
explosion HST observatio ns and find exac tly the same trans- 
formed SN position as in ILi et al.l l|2007h . We also find the 
same point source still visible in the residual image after 
performing PSF-fitting photometry using HSTphot. How- 
ever, we measure the centre of this point source to be some 
~63 mas from the SN position, which given our total as- 
trometric error is a 2.5<r separation. This casts significant 
doubt on the identification of this object as the SN pro- 
genitor. Furtherm ore we a r e una ble to reproduce the pho- 
tometry results of lLi et all (|2007r ) by forcing HSTphot to fit 
at the transformed SN position. Rather we find detections 
of the highest significance (~6.0a in F814W and ~4.4a in 
F450W) when we force a fit at our own measured position 
of this point source, which as we have already said is not 
coincident with the SN position. A more detailed disc ussion 
of this analysis are presented in ICrockett etail l|2009l ). 

Since we cannot confirm that this object is the progen- 
itor of SN 2006ov we derive a 3<r detection limit of mF8i4W 
= 2 4.2. 

ILi et al.l (|2007h derive a mean distance modulus for 
NGC4303 (M61) of fi = 30.5 ± 0.4 (or d = 12.6 ± 2.4 Mpc) 
from two Tully-Fisher distance estimates, and that value is 
adopted here. We apply only a correction for the Galactic 
extinction of E(B — V) = 0.022. Again assuming that the 
progenitor star was a red supergiant we derive a luminosity 
limit of log L/Lq = 4.09 ±0.2 and an 84 per cent confidence 
limit of log L/Lq = 4.29. From Fig. [T]this implies an upper 
mass limit of IOMq . 



The death of massive stars 15 



5.19 2007aa 

NGC4030 is an Sbc spiral with no study of its H II regions 
published. The same arguments as in Section loTTl can be used 
(NGC4030 has Mb = —20.7) to infer an oxygen abundance 
at the galactocentric radius of SN2007aa (0.91r-2s) of ap- 
proximately 8.4 dex. 

The pre-explosion site of SN 2007aa was imaged using 
WFPC2 (2001 July 30) with exposure times of 460s in each 
of the F450W and F814W filters. We determined the posi- 
tion of SN 2007aa on these images using a ground-based 
image of 0.8 arcsec resolution taken with the AUX Port 
camera on the William Herschel Telescope on 2007 March 
11 in the /-band filter. Alignment of this image with the 
pre-explosion frames produced an error circle of 0.07 arc- 
sec for the SN position on the F450W and F814W images. 
No object was detected within this region and hence a 3<r 
detection limit of mFsuw = 24.44 was determined. A kine- 
matic distance estimate of d = 20.5 ± 2.6 Mpc for NGC 4030 
was calculated from its recessional velocity (Virgo infall cor- 
rected) as recorded in LEDA. Extinction du e to the Milky 
Way is estimated to be E(B - V) = 0.026 jSchlegel et al.1 
1998). Assuming the progenitor was a red supergiant we find 
a luminosity limit of logL/Z/Q = 4.38 ± 0.15 and an 84 per 
cent confidence limit of logL/Z/Q = 4.53. From Fig. [1] this 
implies an upper mass limit of 12Mq. 
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Figure 2. The initial mass compared with the final luminosity 
of the STARS and Geneva stellar models. For each mass we plot 
the luminosity at the end of the model, just before core-collapse. 
For the STARS models this is up to the beginning of neon burn- 
ing. The old Geneva models end after core carbon burning. For 
the newer Geneva models both end at silicon burning. The grey 
shaded region represents the range of luminosity for the STARS 
models from the end of core-helium burning to the luminosity at 
the on-set of core neon burning (see Sect. [3] and Fig.fTJ 



5.20 2008bk 

SN2008bk was recently discovered in NGC7793 a nearby 
galax y with two distance modu lus estimates of fi = 27.96 ± 
0.24 (jKarachentsev et all I2T303T ) determined from the tip of 
the red giant branch and /1 = 28.01 from a Tully Fisher 
distance in LEDA. We will adopt the distance modulus of 
H = 27.96 ± 0.24 (d = 3.9 ± 0.4 Mpc). The galaxy is part 
of the oxygen abundance gradient study of iPilvugin et aU 
l|2004| j. The galactocentric radius 0.66r25 and the abundance 
gradient imply an oxygen abundance of between 8.2-8.4 dex 
at the SN position (Mattila et al. 2008), hence we adopt the 
LMC-type metallicity. 

The galaxy has a wealth of prediscovery images avail- 
able from the Very Large Telescope (VLT), with optical BVI 
images from FORS1 and NIR JHK images from HAWK-I 
and ISAAC. We have shown in a recent letter that SN2008bk 
is exactly coincident with a bright, red, source detected in 
the IJHK bands, using high resolution ZG-images from 
the NACO system (Mattila et al. 2008). Although the fore- 
ground extinction toward the galaxy is low and the early ob- 
servations of the SN appear to show no signs of significant 
extinction, Mattila et al. (2008) show that the progenitor 
SED can be fit with a late-type M4I spectral type with a vi- 
sual extinction of Ay = 1.0 ± 0.5. Using methods which are 
entirely consistent with our approach in this paper, Mat- 
tila et al. (2008) have estimated the luminosity and mass 
of this red supergiant. The distance of 3.9T0.4 Mpc and 
A v = 1.0 ±0.5 results in M K = -9.73 ±0.26. Levesque et al. 
(2005) show that using Mk to determine Mboi is preferable 
to using the optical bands. The best fit SED of around M4I 
would correspond to Teff — 3550 ± 50K and a bolometric 
correction BCk = ±2.9 ±0.1 (both from the Levesque et al. 
scale). This results in logL/LQ — 4.63±0.1. and from Fig.[T] 
(at LMC metallicity) this corresponds to a mass of 9t?M Q . 
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Figure 3. Similar to Figure [2] but with initial mass versus effec- 
tive temperature at the end of the stellar model. 

6 SYSTEMATIC UNCERTAINTIES IN THE 
DETERMINATION OF STELLAR MASS 

6.1 Stellar evolution models 

The luminosity estimates and limits used to determine stel- 
lar masses are obviously dependent on the stellar model 
used. Hence it is necessary to compare our stellar models 
to other contemporary models of massive stars. The con- 
stituent physics of modern codes is mostly identical, using 
the same nuclear reaction rates and opacity tables. The dif- 
ferences come from the adopted mass-loss rates, the numeri- 
cal schemes employed to solve the stellar structure equations 
and the tr eatment of mixing, convection and rota tion in the 
codes (see lLanger fc Maederlll995l ; IStanclifre]|2006l. for exam- 
ple) Here we illustrate the differences between our models 
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Figure 4. Similar to Figure[2]but with initial mass versus lifetime 
to the end of the stellar model. 

and those of ISchaller et all <|l992l l. iHirschi et al] (|2004T ) and 
iHeger fc Langerl (|2000l ) . We consider three details of the end 
points of the stellar models, these are the final luminosity, 
the final effective temperature and the stellar lifetimes in 
Figures [2] [3] and H 

One detail to note f irst is that the Geneva rotating mod- 
els (|Hirschi et alj [20041 ) predict a smaller maximum initial 
mass for red supergiant progenitors of 22 Mq rather than 
the 27 Mq from the STARS models, while the non-rotating 
Geneva models (ISchaller et al ] I1992T ) predict a maximum 
initial mass for red supergiant progenitors of 34 Mq. Be- 
yond these masses the codes predict the stars will end as 
H-deficient WR stars (depending on the mass-loss recipe em- 
ployed) . One other noticeable feature is that only the STARS 
models follow the process of second dredge-up and produce 
massive AGB stars at low masses. This is because the other 
models have been stopped before it could occur. Second 
dredge-up is found at similar masses in o ther codes specif- 
ically designed to foll ow this stage (e.g. iPoelarends et all 
l20Qg| ; |Siesj|2006t 120071 ) . 

In Figure [2] the difference in final luminosities between 
models is illustrated. The model sets with the greatest dif- 
ference with our M — logL/L© relation are the non-rotating 
Geneva and Heger & Langer models, while their rotating 
models have reason able agreemen t with the STARS models 
and with the older ISchaller et all (|l992T l models. The main 
reason for the relationships not being exactly similar is be- 
cause of different assumptions of mixing in the stellar models 
and also the tracks end at different points in the stars' evo- 
lution. For example the old Geneva models have lower lumi- 
nosities than the STARS models because they end after core 
carbon burning while the STARS models progress slightly 
further to the beginning of neon burning and we find the lu- 
minosity grows after core carbon burning. Also the Geneva 
9Mq models end after core helium burning and therefore it 
underestimates the final core mass and luminosity. 

The newer Geneva models differ in the treatment of 
mixing and convection in the models which affects the vigour 
of the nuclear burning in the stars and therefore the luminos- 
ity. For example they use a smaller overshooting parameter 
than the older Geneva models as mixing is now also provided 
by the rotation. Thus the rotating models agree with our fi- 



nal luminosities but the non-rotating star luminosities are 
0.3 dex lower. We emphasise that the new non-rotating mod- 
els are artificially pushed to lower luminosities as the mix- 
ing efficiencies (from overshooting) have been significantly 
reduced. Otherwise employing the same mixing parameters 
as previously employed and adding rotational mixing would 
push all the luminosities too high to be consistent with ob- 
served HRDs. 

In general the uncertainty in final luminosity due to the 
assumption of a certain set of stellar models is typically 0.1 
dex between our STARS models and the most up to date 
rotating models. However the issue of how much mixing is 
included and by which mechanism can lead to an uncertainty 
of up to 0.3 dex. This does not pose a major problem to 
our estimates as we are using the luminosity at the end of 
core helium burning to estimate the 84%-confidence upper 
mass limits. In Fig. 2 one can see that the grey area (which 
highlights the region in the STARS code between end of core 
He burning and the end of the model as discussed in Sect.[3| 
brackets nearly all the tracks. 

While the initial mass to final luminosity is uncertain 
the final helium core mass to final luminosity relationship is 
much tighter. This is because the size of the helium core is 
the major factor in determining a red supergiants luminos- 
ity. To remain consistent with progenitor masses measured 
from cluster turn-off ages and with previous studies we have 
determined initial masses (Mi) for our progenitors. Helium 
core masses (MHe-core) can then be estimated from these 
by using the following relation determined from the STARS 
stellar models: 

M Hc -corc = (1.14 ± 0.05)M ; - (0.016 ± 0.002)M? 

-(0.50 ±0.29). 

(1) 

The surface temperatures in Figure[3]show that the final 
pre dicted effective tempera tures are all within 0.05 dex with 
the lHeger fc Langerl l|200Cf ) models being coolest. At higher 
masses the temperature sharply increases as the hydrogen 
envelopes in these cases are low mass (< 0.5Mq) as the star 
is stripped due to mass loss. These temperatures are highly 
sensitive to the boundary conditions in the stellar models as 
well as the opacities used, so it is not easy to simply identify 
the reason for the differences between the models. But the 
uncertainty (±200 K) is well below the uncertainty in the 
surface temperature implied from spectral types of observed 
SN progenitors (typically (±500 K, from the colour-spectral 
type estimates). 

The stellar lifetimes in Figure U also show close agree- 
ment. The most discrepant are the rotating Geneva models. 
Rotation increases the hydrogen burning lifetime consider- 
ably by mixing fresh hydrogen into the core and extending 
the hydrogen burning lifetime of the star. The increase, how- 
ever, is less than 0.1 dex and therefore masses derived from 
lifetimes (i.e. turn-off masses for 2004am and 2004dj) are 
consistent between stellar models. 

Hence we conclude that the use of different stellar mod- 
els are unlikely to have a significant effect on the estimated 
masses and mass limits we have derived, especially if a single 
method is employed and all masses are derived on a homo- 
geneous scale. Furthermore while the initial masses may be 
somewhat dependent on the choice of single star models, the 
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Figure 5. Histogram of the Ay values adopted for the progenitor 
stars (shaded bars) compare d to Ay values of red s upergiants in 
the LMC and the SMC from [Levesque et al. | 1120061 1 (open bars). 
There are 18 SN progenitors plotted here in the shaded bars 
(2004am and 2002hh were excluded as explained in the text), 
and 73 RSGs from the combined SMC and LMC samples. 



final helium core masses that our initial mass correspond to 
should be reliable. 



6.2 Extinction determinations 

It is likely that our largest source of error comes from the 
extinction that we assume is applicable to the line of sight 
toward each SN. This is not likely to be a simple systematic 
effect that would change all the mass estimates and limits 
by a constant. However we need to consider if we are consis- 
tently underestimating the extinction toward the progenitor 
stars and by what magnitude. 

The extinctions which have been estimated for each of 
the progenitor stars come from several methods. All of these 
suffer from their own uncertainties and problems and in gen- 
eral we favour taking the mean of different results. The ratio- 
nale is that no single method is clearly superior to the oth- 
ers and a mean of several, possibly problematic, estimates 
is better than adopting one. The extinctions have been es- 
timated by some of the following techniques: measurements 
of the Nal ISM a bsorp tion lines and calibrating this us- 
ing iTuratto et alj ((2003) ; comparing the early continuum 
slopes to the well observed and reliably modelled 1999em 
and also to unphysically hot black body continua ; fitting 
stellar spectral energy distributions (SED) to the surround- 
ing massive star population within about 10-100 pc; and if 
the SN exploded within an H II region or compact cluster 
then using the value determined from the nebular emission 
lines or the cluster SED. An example of the applications of 
all of these method s applied to SN2001du can be found in 
ISmartt et all (120031 '). In the latest case of SN2008bk there 
is n o accurate ext i nction measurement toward the SN yet 
and iMattila et al.l ((2008) have employed Ay = 1.0 to fit a 
late type M4I to the observed IJHK progenitor colours. 
An extinction of less than this results in a star which is 
intrinsically too red to be compatible with known massive 
red supergiants. We do not revisit the reliability of every 
method applicable to each event as this is dealt with in the 
relevant references cited in the subsections for the events 
above. However we should consider if these methods as a 



whole are applicable and what are the likely sources of er- 
ror. 

The primary concern is our assumption that the extinc- 
tion toward the SN (and surrounding stellar population) is 
directly applicable toward the line of sight to the progeni- 
tor. Two methods probe the intervening line of sight directly 
to the SN. However the early soft X-ray, UV-optical flash 
of the explosion could conceivably have photo-evap orated 
subst a ntial circumstellar d u st clo se to the progenitor i|Dwekl 
Il985t ). IWaxman fc Draind (|2000h have suggested that the 
X-ray and UV afterglow of a GRB could photo-evaporate 
a large cavity surrounding the progenitor star. Scaling the 
GRB energy to the observed flux from recent shock break- 
outs observed for type II-P and Ibc SNe, Botticella et al. 
(in prep) have estimated how much dust could conceivably 
be destroyed in a dense circumstellar envelope. It would ap- 
pear that it is quite possible for such a UV, soft X-ray flash 
to destroy dust masses that could provide several tens of 
magnitudes of extinction in the optical l/-band. 

This is of obvious concern when one considers that 
the observations of luminous red supergiants in the Mag- 
ellanic Clouds, the Galaxy and the L ocal Group are known 
to produce large q uantities of dust |van Loon et al.1 120051 ; 
iMassev et al J 12003 ). A histogram of extinctions toward op- 
tically selected red supergian ts in the LMC and SMC clus- 
ters bv lLevesque et all l|200rj ) suggests that RSGs tend to be 
redder (by on average Ay ~ 0.4 m ) compared to the extinc- 
tions toward the other OB-stars in their stellar associations. 
The mean extinction toward LMC and SMC RSGs is 0.60 m 
and 0.73 m respectively. The mean extinction that has been 
determined toward our SN progenitors is 0.7 ± 1.1 (from Ta- 
ble^ and the large standard deviation is due to 2002hh with 
Ay — 5.2. In this calculation we have left out 2004am which 
clearly suffers from high extinction in line of sight in M82; its 
host cluster is heavily reddened hence its high Ay is unlikely 
to be due to CSM dust shells (see Sect l5.1f|) . This simple 
comparison would suggest there is no clear difference in the 
extinctions of the two samples. If we leave further exclude 
2002hh (as an anomalously high extinction object), we have 
a mean extinction toward the 18 progenitors of 0.44±0.34 In 
Fig. [5] we show the histogram of our Ay estimates towards 
the likely red supergiant progenitors and compare them to 
the LMC and SMC combined population (disregarding the 
two highest values for 2004am and 2002hh). There is some 
evidence to suggest that we have more progenitors in the 
lowest bin — 0.2 m than would be typical for RSG progeni- 
tors. This is not unexpected as for several of our SNe we have 
been forced to adopt the extinction towards the host galaxy 
alone due to lack of additional information. Given the low 
numbers of objects and the differences in the sample size, the 
distribution between 0.2 — 4.0 does not appear to be a major 
cause for concern. There are five events for which we adopt 
a low (foreground Milky Way component only) extinction of 
Ay < 0.1. If these have been underestimated by Ay ~ 0.3 m , 
that would bring the mean Ay of the progenitor sample into 
line with the SMC/LMC RSG populations. In doing so the 
luminosity and mass limits for each event would increase by : 
1999br (logL/L =4.88, < 16M ) ; 2003ie, (log L/L =5.49, 
< 27M ) ; 2006my, (log L/L =4.55, < 13M ) ; 2006ov 
(logL/L =4.35, < 11M S ) 2007aa (logL/L Q =4.6, < 13 
M ). This does not affect the lower mass limit that we derive 
below for the sample from the maximum likelihood analysis 
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and has a minimal affect on the maximum mass as we shall 
see (Sect.0. 

The extinction remains the major source of uncertainty 
and there exist populations of dusty red supergiants which 
are obscured in the visual and near-IR (often by ~ 10 m in 
the V band) and are mid-IR bright as their optically thick 
dust shell is heated by the stellar luminosity and this light is 
preprocessed into thermal mid-IR emission fro m dust grains 
(|van Loon et al.ll2006l. 120051; ILoup et~allll997h. These would 
not appear in the lMassev fc Olsenl (|2003l) and lMassevl (|2002h 
sample as they are too faint optically. However the relative 
numbers of red supergiants (excluding AGB stars, which are 
below the mass threshold to produce SNe) which are visu- 
ally obscur ed (e.g. objects similar to those in lvan Loon et all 
Il998l . l2005h and those which s uffer moderate extinctio ns (the 
optically detectable stars in iMassev fc Olsenl 120031 ') is un- 
known. Such a study to quantify the latest stages in stellar 
evolution in a complete manner would be highly desirable 
and the Magellanic Clouds would appear to be an excellent 
laboratory. Clearly we do see a large popu lation of RSGs 
with l ow- moderate extinctions as shown by iLevesque et al.l 
( 2006) and lMassev fc Olsenl ((2003), and we suggest that our 
progenitors are part of this population. How many dust ob- 
scured RSGs which are missing from optical and near-IR 
surveys remains to be seen. 

Additionally if a mass-loss mechanism (such as pulsa- 
tions) occurs during the final stages of evolution of most 
massive stars as core-collapse approaches one might envisage 
that the progenitors become systematically more obscured. 
This would invalidate the comparison with the LMC RSG 
population. Such severe mass-loss is not well constrained 
observationally or theoretically but if it occurred frequently 
one would expect to see signatures of circumstellar gas as 
well as dust. The type II-P tend to be low- luminosity radio 
and x-ray emitters and tend not to show n arrow hydrogen 
or he lium lines suggestive of CSM shells (| Chevalier et al.l 
I2006D . 

While there is no clear evidence that dense dust shells 
form around II-P progenitors, we at present cannot rule out 
some visual obscuration due to an optically thick dust shell 
which was then evaporated by a soft X-ray, UV and op- 
tical flash at shock breakout. This has been su ggested as 
a possible mechanism for SN2008S j|Prieto et all 120081 ) . see 
Sect.EHl 
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7 MAXIMUM LIKELIHOOD ANALYSIS OF 
THE MASSES OF PROGENITOR STARS 

Using the measurements of progenitor masses in Table[2]it is 
possible to estimate parameters that describe the progenitor 
population. The three parameters of the progenitors that we 
are interested in are the minimum initial mass for a type II-P 
SN, the maximum initial mass and the initial mass function 
(IMF) of the population. 

Estimating these from a small sample is not difficult but 
the relatively small number of data points can restrict the 
accuracy with which one can constrain the most probable 
values. We therefore use the unbinn ed maximum likelihood 
method, e.g. iJegerlehner et al.l 0996). For a large number of 
objects this effectively becomes a \ 2 method. The likelihood 
is defined to be, 



25 - 



£ 20 



15 - 




5 6 7 8 9 10 

Minimum Initial Mass / M 

Figure 6. Plot of the likelihood function for the mass ranges of 
type II-P progenitors. The star indicates the parameters with the 
highest likelihood and the contours the confidence regions. The 
dotted grey lines show the results using the six detections only, 
which results in a lower mass limit of 8.5Mq. The solid black 
lines show the contours using the fixed lower limit and allowing 
the maximum mass to vary. 
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where Pi is the probability of the ith event to have 
mass m. We must define a function for the probability of 
each event, Pi, and then maximize the likelihood to find the 
parameters that give the most probable set of events. 

To make the maximization more straight-forward we 
first take the natural logarithm of the likelihood function 
and so we are required to calculate a sum rather than a 
product, 



ln£ = ^lnP,(r 



(3) 



The probability function that describes the probability 
that a progenitor will have mass m within a certain mass 
range is essentially the IMF. We do need to treat the detec- 
tions and the non-detections differently. For non-detections 
we adopt the probability function, 



Pi = 



(r 



-dm, 



(4) 



where P is the value of the IMF, Salpeter being —1.35, 
yrimin is the minimum mass for a type II-P SN, m max is the 
maximum mass for a type II-P SN, and TO^iimit is the upper 
mass limit for the ith non-detection. If m max is less than 
TOj,iimit we set Pi = 1. If m m i n is greater than m^, limit we set 
Pi = 0.16 because our mass limits are 84 percent confidence 
limits so there is a chance that the progenitor could be more 
massive. This is only important for the mass limit from SN 
2004et. 

For detections the case is more complicated. The er- 
rors for the progenitor luminosity are roughly Gaussian, but 
converting to an initial mass affects the distribution. We 
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first take the best estimated mass, rrij, as the most probable 
value for each detection. Then above this mass we integrate 
the IMF up to the upper uncertainty on the mass estimate, 
'Tij i high- Below rrij we assume the probability distribution is 
a straight line going to zero at the lower uncertainty on the 
mass limit, TOj.iow While these error distributions are some- 
what arbitrary they avoid skewing the overall distribution 
to higher or lower masses as happens when using a Gaussian 
distribution to describe the uncertainties. We have experi- 
mented with different probability functions for detection, for 
example using triangular and rectangular error functions at 
the low and high uncertainty limits. The m max and m m i n 
would typically vary by sa IMq which is within the uncer- 
tainty we derive for these parameters. We feel our chosen 
method for describing the probability function in the case 
of detections is the best representation of the asymmetric 
errors on the mass estimates. Hence for detections we use 
the following probability distribution 



rrij ( m - m j,low) m j 
^j'.low ( m j — m j .low ) ( m mi n 
I r m J,high _ m r ~ 1 

Jm-: 



-dm 



(m 1 



ax) 



-dm. 



(5) 



If mj.iow is lower or higher than m m in and m max then 
the integral is truncated within these limits. 

We calculate the likelihood using the masses for the SNe 
progenitors listed in Table [2] and allow mi ii and m max to 
vary, while fixing the IMF slope to Salpeter (F = —1.35). 
Originally we attempted to let the IMF vary as well as the 
maximum and minimum mass values but find it constrains 
the IMF only very weakly and we chose to fix it at three 
different values, as justified below. 

Furthermore we estimate the confidence regions from, 



ln£„ 



ln£=- X) 



(6) 



where for two parameters when \ — 2.3, 4.6 and 
6.2 we have the 68, 90 and 95 percent confidence regions 
(|Press et al.lll992T ). 

We first estimate m m i n and m max using the six detec- 
tions only (without incorporating the upper limits) . The re- 
sults can be seen in dashed contours in Figure |S] The pa- 
rameters we estimate are m m i„ = 8.5lj 5 Mq and m mal = 
16. o-^O- We then recalculate the likelihoods using both 
the detections and upper limits in the analysis but fix the 
minimum initial mass to 8.5M0 as derived from the detec- 
tions alone. This is because the non-detections only pro- 
vide meaningful information on the maximum initial mass 
m max . In contrast they provide only a weak constraint on 
m m i n as when combined with the IMF they would simply 
favour a low mass due to the rising probability of hav- 
ing more low mass stars. We suggest it is more reasonable 
to calculate m m i n from the detections and effectively this 
converges toward the lowest masses detected in the pro- 
genitor sample. The upper mass limits have a strong im- 
pact on the uncertainty on m max and we determine that 
m max = 16.5 ± 1.5M©. With the error on m max reduced 
significantly, this suggests that at 95% confidence level the 
maximum initial mass to produce a type II-P is 21M©. 

In all of the above we have assumed a Salpeter IMF 
and it is reasonable to question the validity of this assump- 



tion. lElmegreenl (|2008T ) has recently reviewed evidence for 
the variation in the IMF slope in local star forming envi- 
ronments such as Galactic and Local Group Galaxy clusters 
and field populations within the Galaxy and the Magellanic 
Clouds. In clusters and OB associations with total masses 
between 10 2 — 1O 4 M0 there is little evidence for strong and 
real deviations from the Salpeter slope of F = —1.35 above 
and beyond the RMS measurement errors determined in 
each region of high mass stars (with typical uncertainties 
of order ±0.1 to ±0.3). That is not to say that such real 
variations do not exist, only that stochastic affects mean 
that determining the true IMF in localised regions can at 
best reach the accuracy of a few tenths. There is some evi- 
dence that flatter IMF slopes exist in very den se starform- 
ing r egions such as NGC3603 (F = -0.9 ± 0.15; IStolte et al 



2006) and the Galactic centre (r = -1.05 ± 0.05; iKim et al 
20061 ). Also there is evidence that the field population may 
show much steeper slopes, with Y = —1.80 ± 0.09 applica- 
ble for a large sample of fiel d stars lying outsid e clusters 
and associations in the LMC (jParker et al.lll998l). Extreme 
values of around 3-4 have even been found (|Massevl |2002| ; 
iGouliermis et al.ll2002l ) though it is unclear to what extent 
this is simply du e to stellar drift out of low mass clusters. 
lElmegreenl (|2008l ) surmises that F = —1.35 appears to be 
fairly typical in moderate mass clusters and starforming re- 
gions and variation around this is, on the whole, limited 
to approximately ±0.5. As our SNe, and their progenitors, 
seem to reside in typical star forming regions and the field of 
their host spirals there is no compelling evidence to favour 
an IMF too dissimilar to Salpeter. In Fig.[7]we have recalcu- 
lated the ma ximum likel i hood values with the extreme IMFs 
suggested in lElmegreenl (120081 ) of F = -0.85 and Y = -2.0. 
For the shallow IMF slope of Y — —0.85 the best estimates 
of the minimum and maximum initial mass are unchanged 
but the uncertainties increase slightly to m m i n = 8.5^:^0 
and m max = 16.5l 3 M0. This shallow IMF is unlikely to be 
representative of our progenitor environments as they are 
not (apart from perhaps 2004am and 2004dj) in dense clus- 
ters such as seen in NGC3603 and the Arches cluster at the 
Galactic centre. For the steeper IMF, the most likely min- 
imum initial mass increases to m m i„ — and the 
95% confidence limit for m max is pushed to the higher value 
of 22M©. In summary there is no strong evidence (from Lo- 
cal Group studies as reviewed byjElmcgrecn 2008]) that our 
progenitor population should come from a massive stellar 
population with an IMF slope significantly different (i.e. by 
more than ±0.5) than Salpeter, and adoption of the either of 
those extreme values does not significantly affect the values 
of TO m i n and m max . 

Two remaining uncertainties are extinction and the 
value of Hq. As discussed in Sect. 16.21 if we have underesti- 
mated the extinctions toward five events with non-detections 
and replace them with the slightly higher masses, m m i„ and 
m max change by less than O.IM©. As discussed in Sect.0 if 
we employ Ho = 65km s _1 Mpc -1 , then the luminosity dif- 
ferences of the five progenitors for which we employ only a 
kinematic host galaxy distance (see Table[5J) would change 
by ±0.22 dex. This corresponds to approximately 2-3M0 in 
the ZAMS estimate. The value of m m i„ does not change, 
but the the maximum mass increases to m max = 18.5 ± 2. 
This is due to SN1999ev having the most massive progen- 
itor estimate and the host of SN1999ev has a kinematic 
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distance only. Similarly, using Ho = 85km s _1 Mpc -1 keeps 
rrimin within 8-9M0(as all the SNe which determine this 
number have distances from other methods), but the maxi- 
mum mass reduces to m max — 15.5^o 5M0. This illustrates 
that in the future it is important to try to find the type II-P 
SNe from the highest mass progenitors to tie down m ma x as 
reliably as possible. 



8 DISCUSSION 

With our analysis of the progenitor observations and mass 
estimates we are able to consider some outstanding ques- 
tions on the nature of supernovae progenitors from a firm 
observational footing. There has recently been much dis- 
cussion o n the initial masses of progenitor stars of SNe of 
all ty pes (|Gal-Yam et al.ll2007l ; [Li et , al.|[2007l: ISmartt et al ' 



2003) and the nature of fa i nt type II-P SNe JPastorello et al 



2006 ; iNomoto et al. I I2OO3I : IZampieri et alj|2003t Il998) 



Our maximum likelihood analysis reveals that the pro- 
genitors of type II-P arise from stars with initial masses be- 
tween 8.5^1 5 and 16.5 ± 1.5M0. The derivation of the mass 
range assumes that the stars are red supergiants, in that to 
transform the optical or near infra-red limiting magnitudes 
to a luminosity and mass we must assume a stellar progen- 
itor spectrum with a suitable photospheric temperature (or 
range of temperatures). This is well justified in that four 
of the detections have colours consistent with them being 
late K-type to mid M-type supergiants and the requirement 
that a II-P lightcurve results from the explosion of a star 
which has an extended H-rich envelope. (R* ~ 500 — lOOOR© 
lArnettlll980l ; IChevalierl [l976l ; |Popovi ri993). The maximum 
initial mass is important to constrain what the final evolu- 
tionary stage of the most massive stars and the lowest initial 
mass that could produce a type Ib/c, or perhaps II-L and 
Iln, explosion. The minimum initial mass that can support a 
SN explosion is of great interest for explosion models, stellar 



evolution, comparing with massive WD progenitor masses 
and galactic chemical evolution. 

The mass range that we find for the progenitors is 
much low er than ejecta masses of a sample of II-P SNe sug- 
gested by lHamuvl l|2003h . This study estimated ejecta mass 
of bet ween 14-56M0 from the application of the iNadvozhinl 
(2003) formulae to determine energy of explosion, radius of 
progenitor and ejected mass. Even though the error bars on 
the masses are large there is a clear discrepancy between 
our results. The determination of the ejecta masses is very 
sensitive to how the mid-point in the lightcurve is defined to 
determine V50 (the visual magnitude at 50 days) and V50 (the 
ejecta velocity at the same point). The measurement of the 
latter is also highly dependent on which i onic species i s used 
to measure the photospheric velocity and INadvozhinl (|2003l ) 
suggests that the bolometric lightcurve should be used to de- 
fine the plateau mid-point. It appears to us that the choice of 
the point at which to define the measured parameters has a 
critical effect on the physical values determined and caution 
shou ld be employed w hen applying this me thod. We note 
that INadvozhinl <|2003h . with similar data to lHamuvl (|2003f ) 
has determined ejecta masses in the range 1O-3OM0, closer 
to our progenitor mass range but still systematically higher. 
It is likely t hat the t he mid -points of the plateau lightcurves 
defined by lHamuvl l|2003l ) (and the parameters thus aris- 
ing) were not exa ctly compatible with those required for the 
iNadvo zhin (2003) equ ations to be applied. The lower ejecta 
masses of INadvozhinl (120031 ) are probably more reliable in 
that they are estimated with the appropriate input param- 
eters and are a better match to the progenitor masses we 
determine. 

One may ask if a ZAMS mass of 8.5M0 is large enough 
for a long plateau phase to be sustained. In our model the 
star would loose O.5M0 due to stellar winds and with a 
neutron star remnant of 1.5M0, this leaves about 6.5Mo for 



the ejected mass. lHendrv et"aL (2006, 2005) showed that the 
low progenitor masses of SN2004A and SN2003gd (8-9M©) 
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were consistent with the observed recombination powered 
plateau duration, but only just within the error bars of both 
model estimates (see also ISmartt et al.l l2003T j In a future 
paper we will analyse the lightcurves of a large subset of the 
SNe presented here to determine if their progenitor mass 
estimates are consistent with the ejected masses required to 
produce their plateau phases. 

8.1 The minimum mass of II-P progenitors 

Theory predicts that a few of the low-mass progeni- 
tors should be massiv e AGB stars, sometime s referred to 
as Super- AGB stars [lEldridge fe ToutI l2004bl : ISiessI 120071 : 
iPoelarends et all I200S ). The cores of these objects never 
reach high enough temperatures to produce iron, rather the 
oxygen-neon core grows to the Chandrasekhar mass and 
an electron capture SN is triggered. These explosions have 
been predicted to prod uce less luminous SN than in nor- 
mal iron-core collapse (jKitaura et al 1 l2006h . and perhaps 
this signature could be used to find real candidates and 
to identify progenitor stars at the lowest mass range (see 
the Sect . 18.41 and references therein for a discussion of the 
lowest luminosity SNe) . From their models of super- AGB 
stars, iPoelarends et all (|2008T ) suggest that the number of 
these stars at solar metallicity would result in them pro- 
ducing 3 percent of the local core-collapse SNe. This in- 
creases to greater than 10 percent at metallicities below a 
tenth solar. From the observational properties, one of the 
best studied examples of a low-luminosity, low ejecta veloc- 
ity event is SN2005cs, and indeed we do suggest it had a 
low progenitor mass of 8 ± 2M0 . However Eldridge, Mattila 
& Smartt (2007) show that there is a clear observational 
signal for AGB stars, in that these progenitors should be 
much cooler than higher mass red supergiants and hence 
be quite bright at near infra-red (NIR) bands. Deep NIR 
pre-discovery images were available for SN2005cs, and we 
showed that it was unlikely to be a massive AGB star. Thus 
we suggest that all of the 20 progenitors were genuine Fe 
core-collapse events, and we have no evidence for any of 
them being electron-capture events in ONe cores. We also 
have no evidence to support the idea that stars in the range 
~7-9M0 go through 2nd dredge-up and end as quite high lu- 
minosity progenitors, either as S-AGB stars with ONe cores 
or genuine Fe core-collap se events. Our models (and those 
of IPoelarends et al.l 120081 ) would suggest that the luminos- 
ity of these events can reach 4.6 ^ loglz/L©^ 5.1 which is 
significantly more luminous than any of t he progenitors de- 
tected so far. lEldridge fe Tout! (|2004bh and lPoelarends et all 
(2008) point out that this evolutionary phase is very depen- 
dent on the treatment of semiconvective mixing and con- 
vective overshooting. The fact that we don't see luminous 
progenitors with log L/~Lq> 4.6 dex (the highest of our sam- 
ple : SN2005cs) would apparently disfavour the scenario in 
which ~7-9M0 stars increase their luminosity due to 2nd 
dredge-up before collapse. The apparent luminosity of the 
progenitors (around 4.3 — 4.6) favours a lower limit than is 
normally assumed for core-collapse with no luminosity spike. 

The lower mass limit we derive from Sect .[7] of 
8.5~^\ 5M0 is interesting to compare with the maximum stel- 
lar mass that produces white dwarfs. A compilatio n of mass 
estimates of white dwarfs by iDobbie et all (2006) suggests 
that, in Milky Way intermediate age clusters, stars up to 



6.8-8.6M0 produce white dwarfs and the y suggest this as the 
initial mass range for core-collapse SNe. lRubin et al.l (|2008T l 
suggest that a homogeneous analysis of WDs in their Lick- 
Arizona White Dwarf Survey (LAWDS) confidently deter- 
mines the maximum mass to be no less than 6M0. Further 
recent e vidence suggests tha t the mass limit is no less than 
7.1M (| Williams et a]||2009h . A slightly higher mass limit is 
not ruled out as there is ongoing work on younger clusters to 
find WDs and determine their progenitor age (C. Williams, 
private communication). Hence the two, very different ap- 
proaches, of SN progenitor mass and WD progenitor masses 
appear to be converging toward 8± lM© . Unless both meth- 
ods are significantly in error it would seem unlikely that the 
lower mass for a core-collapse SN is outside this range. The- 
oretically several ma ss-limits have been determined, rang- 
ing from 6 to HMm dRitossa et al.lll999l; [Heger et alj|2003l ; 
lEldridge fe ToutI l2004bl ; IPoelarends et al .1120081 . and refer- 
ences therein) depending critically on the amount of con- 
vective overshooting employed. In our analysis we have used 
models with convective overshooting as there is growing evi- 
dence that extra mixing is required above that predicted by 
mixin g-length theory fe.g. ISchroder et al.lll997l ; lAerts et al.l 
2003). We suggest that a minimum initial mass of lOM© or 
more can now be ruled out for two reasons. Firstly we detect 
four progenitors with best estimated masses below lOM© 
although admittedly the individual uncertainties would not 
rule out a higher mass progenitor. Secondly in our maximum 
likelihood analysis masses at lOM© and above are ruled out 
at over 95 percent confidence, even with a steep IMF of 
r = —2. This value is supported by the fact that type II-P 
SNe are not always associated wi th underlying H 11 emission 
line regions in their host galaxies I Anderson fe Jamesl (|2008T l 
which would suggest their progenitors are from a population 
of less than ~10M©. 

We suggest that 8.5 jl^ 5 Mq is the current best estimate, 
based on observational constraints, for the lower limit to 
produce an Fe core-collapse driven SN of type II-P. This 
is in agreement with the mass range for the most massive 
progenitors of WDs. 

8.2 The maximum mass of II-P progenitors 

The maximum mass of a star that can produce a II-P su- 
pernova is an important threshold to constrain. In Fig. 
[5] we summarise the initial masse s of all the p rogen itors- 
Similar plots we r e first shown by iHeger et al.l (I2003T ) and 
lEldridge fe ToutI l|2004bl ) with a large range of metallic- 
ity plotted on the vertical axis, from super-solar to metal 
free. As our progenitor stars cover a relatively small range 
in metallicity, we have removed the axis scale and in- 
stead flagged the points with a metallicity coded grey scale. 
Clearly the highest mass of a detected progenitor is W^M© 
(SN1999ev) with one upper limit above 2OM0, due to shal- 
low pre-explosion images (SN2003ie). Our estimated maxi- 
mum initial mass for a II-P (Sect.[7)l is m max = 16.5±1.5M0 , 
with a 95% confidence limit (assuming Salpeter IMF V — 
— 1.35) of 21M©. Fig.[5]is effectively a cumulative frequency 
distribution (CDF) which is constrained at the lower and 
upper mass limits and has an IMF with V — —1.35 consis- 
tent with the limits in between (the CDF of the Salpeter 
IMF is plotted as the thick grey line). This Salpeter IMF is 
a good fit to the distribution of masses and mass limits, if 



22 S.J. Smartt et al. 



1 999ev 
2004dj 
2003ie 
2002hh 
1999an 
2004am 
1999em 
2001du 
1999br 
1999gi 
2007oo 
2006my 
2004dg 
2006bc 
2006ov 
2004et 
2008bk 
2005cs 
2003gd 
2004A 



- ■ 



i — ■ 



-ZH 
-H- 



15 20 
Initial moss / 



Figure 8. The initial masses of all our type II-P progenitor stars, 
compared with our theoretical limits for production of supernovae 
of different types and type of compact remnant. The box symbols 
are shaded on a metallicity scale, the lighter the shade the lower 
the metallicity, with the values taken from Table[2] 



the hard minimum and maximum masses for II-P progeni- 
tors hold. Stars more massive than about 2OM0 would be 
easily detectable in our archive images, and there is unlikely 
to be any bias against detecting the most massive progen- 
itors. Hence there does appear to be a real upper limit to 
the mass of stars that produce normal type II-P SNe. The 
one caveat to this is if the progenitor stars suffer large cir- 
cumstellar extinctions which are photo-evaporated in the 
explosion. We discussed this in Sect.[6]and while we cannot 
see a compelling case for such an effect in our population we 
cannot rule it out. 

We can compare this maximum mass limit with the ra- 
tios of CCSN types in Table[ T1 Wit h a maximum possible 
stellar mass of 150M© (Figer 2005), the fraction of stars 
born with masses between 8.5-16.5M0 (for a Salpeter IMF, 
T = —1.35) is ~ 60 per cent, closely mirroring the type II-P 
rate. One might immediately conclude that the agreement 
suggests that all stars above ~17M0 produce the other va- 
rieties of CCSNe. However this is too simplistic and ignores 
our wealth of knowledge of massive stellar populations from 
Local Group studies and interacting binaries. 



8.2.1 The red supergiant problem 

Massive red supergiants have been frequently surveyed in 
the Milky Way and the Magellanic Clouds, and up un- 
til recently their luminosities as determined from model 
atmospheres implied th at they are found at evolutionary 
mass es up to 4O-6OM (|Massev fc Olsenll2003l ; iHumphrevsl 
19781). Howeve r usin g new MARCS atmosphere models 



Levesque et alj (|2006l ) have shown that the effective tem- 
peratures of these stars have been revised upwards and they 
have combined this with revised bolometric luminosities 
based on A"— band magnitudes. The result is that the highest 



luminosity red superg iants of iMassev fc Qlse3 (2003) and 
iLevesque et alj (|2005h now have warmer effective temper- 
atures and luminosit i es tha t imp ly masses of bet ween 12- 
3OM0. IMassev et all (|200lh and ICrowtherl l|2007l) suggest 
stars with an initial mass of around 25M© could evolve to 
the WN phase in Galactic clusters, at solar metallicity. The 
mass estimates generally come from the estimated age of the 
stellar clusters as measured from the turn-off. Only two out 
of 11 in the Massey et al. study are as low as 20-25Mq and 
one can really only take this as a lower limit. The minimum 
initial mass to form a WR star in the LMC (and SMC) 
has been estimated at 30Mp) (and 45M q respectively) us- 
ing similar methods l|Massev et al,ll2000l ). Stars above these 
masses, if they explode as bright SNe, should produce H- 
deficient (and He deficient) SNe like the Ib/c we observe. 
Hence there is good agreement between the maximum ob- 
served masses of red supergiants in the Galaxy and the LMC 
and the minimum mass required to produce a WR star, from 
the ag es and turn-off masses of coeval clusters. ICrowtherl 
(2007) points out that there are few Milky Way clusters 
that harbour both RSGs and WR stars which would sug- 
gest that there is a definite mass segregation between the 
two populations. The metallicity ranges of our progenitor 
sample (Table[2} range between solar and LMC, hence these 
studies of Local Group stellar populations would suggest 
the minimum initial mass for a single star to become a WR 
(probably of type WN) is 25-30M o . 

The question is what is the fate of the massive red super- 
giants between 17Mq and 25-30Mq ? They appear to exist 
in this mass range and one would expect them to produce 
SN of type II-P but they are missing from our progenitor 
population. A single star of initial mass of 17M© does not 
have a high enough mass-loss rate to strip its outer layers 
of enough mass to become a WR star and hence a lb or Ic 
SN (either observationally or theoretically). If our sample 
of 20 progenitor stars were really sampled from an underly- 
ing population of red supergiants, with initial masses in the 
region 8.5 — 25M©, then a Salpeter IMF would suggest we 
should have 4 between 17 — 25M©. The probability that we 
detect none by chance is 0.018 (or 2.4cr significance). For a 
steeper IMF of V = —2.0 the numbers are 3 stars, proba- 
bility of 0.05 and 2a significance. We term this discrepancy 
the "red supergiant problem" , in that we have a population 
of massive stars with no obvious channel of explosion. 

One could attempt to fill this mass-gap with the other 
SN types Iln and II-L and lib. The fraction of stars born 
with masses between 17 — 25M© (within an underlying pop- 
ulation of 8.5 — I5OM0) is 18 per cent, and Table[T] suggests 
the combined rate of II-L, Iln and lib is 12 per cent. Hence it 
is perhaps appealing to account for the red supergiant prob- 
lem by saying that at least some of these stars form II-L, 
Iln or lib S Ne. B ut there is evidence arguing against this. 
iThompsonl (|l982l ) presented a deep photographic plate of 
NGC6946 49 days before the maximum of the II-L SN1980K 
and found no progenitor or discernible stellar cluster. He 
suggested an upper mass limit of < 18Mq and using our stel- 
lar tracks and more rec e nt dis tance we recalculated this as < 
2OM in lSmartt et all (120031 ). SNe lib have been suggested 
to be from interacting binary systems and for SN1993J a 
viable model is a close pair of 14 and 15M© stars. The bi- 
nary companion to SN1993J's progeni tor was theoretically 
predicted and observationally detected (jPodsiaklowski et al.l 
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ll993l:IWooslev et alj[l99i ; iMaund et"al]|2004 ). [Ryder et all 
(2006) suggest a similar scenario explains their detection of 
a stellar source at the position of the lib SN 2001ig. A single 
28M0 WNL star was f avoured as the progen itor of the re- 
cent SN lib 2008ax bv lCrockett et all ([20081 ) . but a binary 
system cannot yet be ruled out. There is also evidence that 
very luminous type Iln arise from very massive LBV type 
stars, generally thought to be > 4OM0 and hence too high 
mass to solve the problem (see Sect 18.2.%) . 

One could appeal to rotation as a way o ut and the 
rotating Geneva models of iHirschi et all i[2004l) predict an 
upper mass limit of 22M© for a hydrogen rich progeni- 
tor (for stars rotating initially at 300 kms -1 ). Above this 
mass a single star ends its life as a WR and hence a Ib/c 
SN. This is well above our estimated maximum initial mass 
17Mq, but consistent with the 95% confidence limit. How- 
ever this would mean every II-P progenitor would have to 
be rotating initially with speeds around 300 km s . This 
is clearly not what we see in the rotational v elocity distri- 
butions in the Galaxy, or Magellanic Clouds. dHunter et ahl 
120081 ; iDufton et al.ll2006l ; iHuang fc Giesl |2006| ), which sug- 
gest less than 5 per cent of massive stars should be rotating 
at such intrinsic rotational velocities. 

8.2.2 Black hole formation 

An intriguing possibility is that the red supergiant prob- 
lem is due to the vast majority of high-mass stars above 
17M© collapsing to form black holes and either very faint 
supernovae or no explosion at all. Theoretically this has 
be en suggested for some time, fo r exa mple most recentl y 
bv lFrverl |l999l ); lFrver et all (|2007l ) and lHeger etafl <[2003l ). 
Our model stars in the mass range of 20-27Mq end as hydro- 
gen rich red supergiants with helium core masses of > 8M0 
and such masses have been suggested to result in the forma- 
tion black holes (this line i s plotted for refere n ce in Fig. [5] 
lFrverlll999l ). The models of iLimongi fc Chieffll (|2007l . 120031 ) 
suggest the maximum mass to produce a type II-P SN is 
3O-35M0 and a minimum initial mass for black- hole forma- 
tion is 25-30M Q . Although fFrveri [|l999b notes that the mass 
range to produce black holes is theoretically quite uncer- 
tain. For example reducing the mean neutrino energy by 
20% could reduce the explosion energy by a factor 2 and 
push the minimum mass for black hole formation to as low 
as 15M . 

As pointed out by IKochanek et all l|2008h . the collap- 
sar model in which a GRB is produced along with a type Ic 
SN, is likely to be too rar e to produce the bulk of the black 
holes seen in our Galaxy [[MacFadven fe Wooslevl[l999h . Al- 
though the collapsar scenario would have problems within 
massive hydrogen rich progenitors (the j et would have dif- 
ficult y in escaping f rom a red supergiant iMacFadven et ahl 
120011 ). lYoung et all (|2005l ) suggest that bright II-L SNe may 
be black hole forming events, in which the collapsar mecha- 
nism occurs within a massive H-rich star. 

Whatever the explanation we have evidence for a lack 
of progenitors above 17M0 and perhaps the minimum mass 
to form a black hole could be as low as this. It could be that 
stars in the 17-3OM0 range produce SNe so faint that they 
have never been detected by any survey. In this case they 
would typically be fainter than Mr ~ —12. If the limit for 
black hole formation is low then it bodes well for surveys for 



disappearing stars. IKochanek et al.l (|2008l ) have suggested 
that a survey of nearby galaxies over several years would 
have a chance of detecting massive stars that disappear 
without an accompanying SN. From a similar comparison 
of a Salpeter IMF with the general progenitor compilation 
of lLi et all (|2007l > they also suggest there may be a dearth 
of massive star progenitors. Their calculation is somewhat 
inexact in that it includes 1999gi and 2001du as possible 
detections and is neither volume or time limited to min- 
imise biases on SN and progenitor selection effects, and the 
masses come from many inhomogeneous methods. But it 
does support our quantitative mass range estimate for II-P 
progenitors. 



8.2.3 Binaries and Ibc SNe 

A further flaw in the argument that the type II-P rates 
match the mass range of 8.5-16 .5M0 is that it ignores the 
consequence of binary evolution. IPodsiadlowski et al.l (|l992l ) 
suggested that around 15% of SNe could be from interacting 
binaries in which mass transfer causes the primary to loose 
its H (and He) envelope. This assumes that about 30% of all 
massive stars are in close binaries that will interact in case A, 
B or C mass transfer. Recent results s uggest that ^ 60% of 
massiv e stars c ould be in c l ose bi naries iKobulnickv fc Frverl 
(|2007h leading [Fryer et"al] l|2007h to claim that perhaps all 
local Ibc SNe could be formed in binary systems and the pro- 
genitors could thus have initial masses down to our m m i n 
limit of 8.5M0. It appears very likely that at least some 
Ibc SNe are formed i n moderate mass int eracting binaries 
l|Crockett et al.ll2007i ; lFddridge et al J 120081) 

At aro u nd so lar metallicity iFrver et alj (|2007l ) and 
iHeger et al.l l|2003l ) argue that single, massive WR stars all 
have core masses large enough to form black holes and 
that they can't be the progenitors of the local, normal 
Ibc SN population. They suggest that these should give 
weak SNe or no explosion at all. Current observations have 
not yet confirmed that massive WR stars are definitely 
the progenitors of Ibc SNe. The bulk of the population 
may form black holes with no explosion and a fraction 
(with low metallicity and high rotational velocities) may 
form black ho les in the collapsar mo del with an accom- 
panying GRB ([Wooslev fc Bloomll2006l ). The ejecta masses 
and pre-explosion limits of SNe Ibc (with no associated 
GRB) are consistent with them being stars of ~1O-2OM0 
strip ped of their envelopes through close b i nary interac- 
tion dMazzali et alj|200a ICrockett et a.1.1120071 : IValenti et ail 
2008). We don't yet have a firm confirmation of a Ibc ex- 
plosion (with no associated GRB) directly associated with 
a massive WR star, or with ejecta masses high enough to 
suggest a WR progenitor. Those broad lined SNe with a 
GRB associated do have high enou gh ejecta masse s to be 
consistent with LMC type WC stars (|Crowthenl20071 ). How- 
ever there is a suggestio n of a 28M0 WNL progenitor for 
the type lib SN2008ax (|Crockett et al.ll2008l ) and this SN 
appears to be o n the H abundance c ontinuum between lb 
and lib events (|Pastorello et al.ll2008l ). We will discuss the 
Ibc progenitor scenarios further in the second paper in this 
series. 
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8.2-4 The type Iln population and their progenitors 

There is clear evidence now that some very massive stars, 
above 25-30M© do explode as very bright SN, SN2005gl is a 
type Iln at 65 Mpc and has a very luminous progenitor de- 
tected at M v ~ -10.3 l|Gal-Yam et al.l 120071 ). The latter is 
evidence that luminous blue variables are direct progenitors 
of some SNe and this is supported by s tudies of the energies 
and spectral evolution of other events ( Kotak fc Vinkl 20061 ; 
ISmith et al.ll2007l : iTrundle et aI1l200Sl : ISmith et al.ll2008bh . 
The case of SN2006jc showed that an LBV-like outburst oc- 
curred directl y coincident with a pe culiar type of hydrogen 
deficient SN fPastorello et al.ll2007f h SN2006jc resembles a 
type Ic with narr ow lines of He ar i sing from a circumste l- 
lar He-rich shell (iFolev et all 120071 ; IPastorello et al.l l2008h . 
IWooslev et al.l l)2007l ) have sug gested that both t he super- 
bright Iln events (2006gy like; (Smith et al.1 120071') and the 
double outburst events (2006ic-like; IPastorello et "all 120081 ) 
may not be the canonical core-collapse mechanism, but be 
due to pulsational pair-instability in the cores. 

Thus one might venture that above 17M0 the vast ma- 
jority of stars form black holes at core-collapse and can't 
produce bright explosions through the canonical neutrino 
driven convection mechanism. A fraction of them form col- 
lapsars due to a co mbination of rotation or binarity and 
low metallicity (see IWooslev fc Blooml [20061 ). And a frac- 
tion may form H-rich luminous type Iln SNe through the 
pulsational pair-instability mechanism. 

A caveat to t his is the discovery of neutron stars i n 
two young clusters ijMuno et al.ll200d : lMessineo et al.ll2008f ). 
which suggests the progenitors had initial masses greater 
than 40Mq and 2O-3OM0 respectively. These stars should 
perha ps have formed black holes but iBelczvnski fc Taaml 
(2008) suggests that under certain conditions, binary evo- 
lution could result in stars as massive as 5O-8OM0 ending 
up as neutron stars. A further argument against is that the 
locations of Ibc SNe tend to be more closely associated with 
Hn regions than II-P SNe, sugg esting a higher initial ma ss 
range for the progenitors of Ibc (jAnderson fc Jamea 2008) 

The nature of the deaths of the most massive stars, 
whether in black-hole forming events, or other explosion 
mechanisms, still remains to be determined. The combina- 
tion of studies of direct progenitor detections, environment 
evaluation, SN ejecta and remnant properties will be a fertile 
field for discovery for many years to come. 



8.2.5 The progenitor of SN 1987a : Sanduleak -69°202 

Although the progenitor of SN 1987A is often quoted to be 
a 20Mq star, one needs to be careful with a simple interpre- 
tation of placing the progenitor on an HRD and taking the 
close st mass track. The spe ctral type and UBV magnitudes 
from Walb orn et al. (1989) suggest a B3-type supergiant 
(T e ff ~ 15750, from t he calibration of LMC B-supergiants 
in ITrundle et al.1 l2007l ) and hence log(L/L ) = 5.1 ±0.1. 
When placed on single-star evolutionary tracks this lies close 
to a 2OM0 model just after the end of core-H burning. 
However it is not valid to assume 2OM0 as the progeni- 
tor initial mass, as the model track is not at its endpoint 
and is no where near to having an Fe-core (or at least at 
the point of neon burning within a helium core). The lu- 
minosity of the He-core of an evolved massive star deter- 



mines the stars luminosity and we estimate the correspond- 
ing He-core mass to be 5.2^i o-^©- The initial mass of a 
single star required to produce t his core mass is 
The interacting binary m odel in iMaund et all |2004) and 
iPodsiaklowski et al.1 i| 19931) can produce a SN1987A like pro- 
genitor with a pair of 14 and 15M0 stars. And a merger in- 
volving a lower mass star of 3-6M0 with an ev olved 15-16M0 
prim a ry can also account for the lum inosity (jPodsiadlowskil 
1 19921 ; iMorris fc Podsiadlowskl 120071 ). In both scenarios an 
initially less massive star gains mass to explode with a final 
mass of 2OM . However the helium core mass was that ex- 
pected for a 15M0 star leading to its position in the blue 
part of the HR D. The ejecta mass has also been estimated at 
around 15M Q l|Arnettlll996bl ). Hence our suggestion of black 
holes coming from ~17M0 stars and above is not directly 
disproven by the example of Sanduleak —69° 202 



8.3 Explosion mechanisms and production of 56 Ni 

The tail phase of type II-P SNe are thought to be powered 
by the radioactive decay of 56 Ni and recent studies have 
shown that there can be a large range in tail phase lumi- 
nosities. This would imply that a different mass of 56 Ni has 
been ejected in the explosions. As 56 Ni is created by explo- 
sive burning of Si and O as the shock wave destroys the 
star, it can b e used as a probe of t he explosion mec h anism . 
For example iTuratto et all (|l998l ). IPastorello et al.l i|2004l ). 
and lNomoto et al.1 (|20o5 T predict that high mass stars may 
undergo fall back in which some of the 56 Ni falls back onto 
a proto-neutron star or black hole and, hence, one might 
get a fainter supernova. This has lead to suggestions that 
plotting initial mass versus ejected mass of 56 Ni could lead 
to a bimodal population. There are estimates of the mass 
of 56 Ni for nine of the SNe in our sample (by ourselves and 
also from other groups, already published in the literature) 
and we can investigate this relation in a direct way. 

The 56 Ni mass can be estimated from the tail phase 
magnitudes using several different methods. For example , 
the bolometric luminosity of the tail phase (Hamuy 2003), 
a direct comparison with SN 1987A dTuratto et alj 1 199 



and t he "steepness of decline" correlation ( Elmhamdi et all 
20031). For t he two re cent SNe 2004A and 2003gd, 
Hendry et al.l l|2005l . 120061 ) have compared the three meth- 
ods and find the first two in good agreement, while the 
mass from the "steepness of decline" relation gives some- 
what lower results (at least for these two events). It is im- 
portant that a consistent method is used to determine all 
the 56 Ni masses if any comparison is to meaningful, partic- 
ularly as the uncertainties on the estimates are often fairly 
large. Hence we determine 56 Ni masses from one consistent 
approach. The values fo r SNe 2004A 2003gd, 1999gi and 
1999em were taken from iHendrvl l|2006l ) who used the bolo- 
metric tail-phase luminosity method of (|Hamuvll2003l ), and 
the distances and reddening alread y adopted in T ableful The 
mass for SN2004dj was taken from IZhang et ail (|2006l ) who 
found that the bolometric luminosity of the tail-phase gave 
a value very similar to that from the "steepness of decline" 
relation. We take their value from the bolometric tail-phase 
luminosity to ensure consistency (0.025 ±O.O1OM0). This is 
very similar to the value determined bv lVinko et al.l (|2006l ) 
with a simple radioactive decay model applied to the bolo- 
metric tail luminosity (O.O2±O.O1OM0). By a similar analy- 
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sis, th e value of SN2004et was determined with the iHamuvl 
i|2003l ) formula, to give 0.06 ± 0.02M Q , the highest of our 
estimates. 

The value for SN2004dg has been determined from the 
late-time, tail phase magnitudes from our HST imaging (see 
Sect l5.12"j) . We determine V = 20.8 and estimate an explo- 
sion epoch from spectra and photometry during the plateau 
phase. Applying the same bolometric tail phase method as 
above results in a value of 0.010 ± 0.005Mpi . A very low 56 Ni 
mass for SN1999br has been reported in faamuvil2003h . al- 
though he used a distance of 10.8 ± 2.4 Mpc, and to keep 
the analysis consistent we scaled his value to our larger 
adopted distance listed in Table[2] resulting in a value of 
0, 003=1=0. OOlM g . This is is in go o d agr eement with the value 
determined by IPastorello et all l|2004h of ~0.002 who used 
a similar distance to ours to calculate the mass from a di- 
rect comparison with the tail phase luminosity of SN1987A. 
The estimates for SN2006ov and SN2006my have been cal- 
culated using a similar bolometric tail phase method as de- 
scribed above to give 0.003 ± 0.002M Q and 0.03 ± 0.015M o 
respectively (Maguire et al., Spiro et al. in prep). 

Finally we consider the case of SN2005c s. The early time 
spectr a have been studi e d exte nsivel y by Pastorello et alj 
l|2006h . iTakats fc Vinkl l|2006h and iBrown et al] (j2007h . 
who all find it to be a moderately faint II-P with low ejecta 
velocities. It appear s a similar type of eve nt to the faint SNe 
1999br and 2002gd (jPastorello et al]|2006l ) , and a measure of 
its luminosity in the tail phase is especially interesting par- 
ticularly as it is one of the events wit h a detected pro g enitor 
and well determined mass (Table©. iTsvetkov etafl ([20061 ) 
present several photometric measurements of SN2005cs in 
the tail phase and suggest a 56 Ni mass of 0.018Mq, which 
would be similar to 2003gd. As the progenitors are likely to 
have been red supergiants of quite similar masses one might 
be en couraged by this agreement. However (|Pastorello et al] 
2009) presents new measurements of the the tail phase 
magnitudes of S N2005cs finds it t o be significantly fainter 
than reported in ITsvetkov et al] (|2006l ). The difference is 
likely due to the differing methods employed to determine 
the luminosity in this faint phase. Pastorello et al. have 
used image subtraction to remove contamination of the host 
galaxy, which can be significant as the SNe fades. These 
fainter measured magnitudes suggest an ejected 56 Ni mass 
of 0.003 ± O.OOIM© , and we believe this to be a more realis- 
tic estimate. For re ference, the det ermined value of O.O75M0 
for SN1987A from lArnettl (|l996al ) is also plotted in Fig. [9] 



8.4 The nature of faint II-P SNe 

In Figure [9] we see how the nickel mass created in II-P SNe 
compares to the initial mass of the progenitor star. Simi- 
lar diagrams have been produced before but these derive 
the initial mass by a model dependent two step process. 
First of all the ejected mass is estimated from modelling 
the SN lightcurve and velocity evolution of the ejecta, and 
then an initial main-sequence mass is inferred by assuming 
an estimated remnant mass and acco unting for the effects 
of mass-loss during stellar evolution (|Nomoto et al] 120031 ; 
IZampieri et al] [20031 ). Such models have given satisfactory 
fits to the luminosity and velocity measurements, but our di- 
rect pre-explosion measurements can provide valuable inde- 
pendent information. It is certainly clear that there is a pop- 
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Figure 9. Plots of initial mass vs. mass of 56 Ni. The grey line 
is the Mq/M (COcore) normalised to pass from the 1987A nickel 
mass. 



ulation of low luminosity II-P SNe, which have lower ejecta 
velocities throughout the photospheric stage and very low 
tail phase bolo metric luminosi t ies. The first one recognised 
was SN1997D dTuratto et all 1 19981 ; IZampieri et~aH Il998l ; 
iBenetti et aL 200 if) and in our sample 1999br and 2005cs 
are similar I Pastorello et alj |2004. 2006]). The clear implica- 
tion of the low tail phase luminosity is that a low mass of 
56 Ni is produced in the explosion. Hence one might hope to 
relate this to the explosion mechanism. 

Two alternative suggestions for the faint II-P SNe have 
been proposed. One is that the SNe formed black holes 
at core-collapse and the 56 Ni produc ed fell back into the 
black hole rather th an being ejected (jZampieri et al]|2003l ; 
iNomoto et al.| [2003). An alternative is that they are intrin- 
sically less energetic explosions of lower mass stars. Stars of 
9-11Mq can have large density gradients in the O-Si lay- 
ers around the proto-neutron star and the shock may pro- 
duce lower temperatures than in higher mass counterparts. 
Statistical equilibrium is only reached in a thinner shell of 
Q and Si rich gas, hen c e low amounts of 56 Ni is produced 
<|Mavle fc Wilsonll 19881 ). fChugai fc Utrobinl (|2000h have also 
favoured this low-mass star scenario in their fitting of the 
nebular spectra of SN1997D wi th a hydrodynamic m odel. 
The high mas s stellar origin of IZampieri et alj (120031 ) and 
INomoto et ail (|2003l ) predicts hydrogen rich progenitors of 
greater than 25Mq, but as shown above we do not detect 
any red supergiant stars of such high masses. These types 
of stars, would be the easiest to detect and in particular 
for 1999br and 2005cs we favour the low mass scenario. 
Figure|9] argues against the high mass scenario for the low 
56 Ni mass SNe, and we find no evidence of the branching of 
the figure at high masse s into low and high energy SNe as 
sugges ted by both lZampieri et alj l|2003l ) and INomoto et all 
(2003) The low 56 Ni SNe have initial masses well below the 
limit required to produce a massive enough core for black 
hole formation. In fact rather than being two separate pop- 
ulations there appears to be a continuous trend with lower 
initial mass pr oducing lower masses of nickel in core-collapse 
l|Maundll2005l ). 

One could advance an argument that the progenitors 
were actually of higher mass and not detected in our im- 
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ages for some reason. The most obvious reason is high dust 
extinction as we discussed in Sect . 16.21 We argued that it 
is unlikely that we have systematically underestimated the 
masses due to large extinctions. In addition, the fact that 
there was no detection of the 2005cs progenitor in deep NIR 
JHK bands argues against a mass ive progenitor surro unded 
in dust, of anything up to Av ~ 5 |Maund et al.| [2005). Also 
SN2008bk was detected in the the NIR {JHK bands) and 
hence even a visual extinction of Av ~ 10 would have only a 
~ l m affect in K, increasing the logi/L© by 0.25 dex. An- 
other possibility is that the stars were not red supergiants, 
but hotter, blue stars. One could invoke this explain the 
non-detections and also argue that the detected red stars 
are heavily reddened bluer objects. But in this case they 
would be more compact and the lightcurves would then tend 
to resemble SN1987A, and there is no evidence such pecu- 
liarities in any of the SNe presented here. If all, or most, 
II-P progenitors are not red supergiants this would pose se- 
rious problems for our stellar evolution theory and models 
of II-P lightcurves. We consider t his possibility unlikely. A s 
this paper was in the review stage, lutrobin fc Chueail (|2008h 
proposed that the ejecta mass of SN2005cs, from hydrody- 
namic modelling of the lightcurve and velocity evolution, 
was around 17.3 ± 1 Mq . They also suggest a progenitor 
mass of 2O-25M0 for SN1999em which is again significantly 
higher than our stellar evolutionary mass. This severe dis- 
crepancy between two methods is interesting and should be 
explored further in the future. 

As the 56 Ni is produced by e xplosive silicon and o xygen 
burning of the mantle material l)Wooslev et alJ 12002). one 
might imagine that the amount of nickel formed in a SN 
is related to the amount of oxygen and silicon in the pro- 
genitors core. In Figure [§] we over-plot the mass of oxygen 
in the carbon-oxygen (CO) core divided by the mass of the 
CO core from our stellar models against the initial stellar 
mass. The line is scaled to fit SN 1987A. This simple model 
follows the general trend of the observations, but with some 
large scatter. Other factors will effect the amount of nickel 
produced in a SN but the relation could suggest that the 
amount of source material is a primary factor in determin- 
ing how much nickel is produced. Alternatively it could be 
a reflection of the density and temperature in the explosive 
burning region in the stars. Whatever the physical reason, 
we have not detected high mass stars as progenitors of the 
faint II-P SNe. Hence we favour the lower mass star progen- 
itor as the origin of these events. 



9 LESSONS LEARNED AND FUTURE 
POSSIBILITIES 

This decade long effort to systematically search and detect 
SN progenitors provides an opportunity to reflect on lessons 
learned and ways to increase discovery potential. It is per- 
haps somewhat surprising that unambiguous detections in 
multi-band images have been difficult. As discussed above 
this could imply that the high mass, H-rich, bright super- 
giants up to around 30Mq are not the progenitors of type 
II-P or II-L SNe. Whatever the reason, it seems clear that 
pursuing this project needs a different approach to make 
firm detections of progenitors rather than providing many 
more upper limits. The latter can be interesting if enough 



are available, but ultimately it is discovery and characteri- 
sation of objects that will advance the field. 

As discussed in Sect. 5, the three progenitors which are 
reliably and unambiguously detected (2003gd, 2005cs and 
2008bk) and have high-quality multi-colour images (giving 
colour information) are all closer than 10 Mpc. Although the 
global image archives of nearby galaxies have been steadily 
increasing over the last 10 years, the total number of galax- 
ies and the quality of the images will not rapidly improve 
in the near future, e.g. it has taken 15 years of post-SMI 
HST operations to get to the current archive content. It 
seems clear that the current HST archive does not contain 
deep enough images of galaxies to enable routine detection 
of progenitors beyond about 10 Mpc. This does not mean 
we should give up on searching for progenitors of SNe be- 
tween ~ 10-25 Mpc as one may be fortunate to detect higher 
mass , and brighter, proge nitors than we so far have done 
(e.g. iGal-Yam et"al]|2007h . Each event should certainly be 
scrutinized for the potential of fortuitous progenitor dis- 
covery but when considering how to actively improve the 
possibilities the following should be noted. In Cycle 10 of 
HST operations (2001 June 01 - 2002 May 31) we were al- 
located a SNAP proposal (SNAP9042) to enhance the im- 
age archive of starforming galaxies within about 20 Mpc and 
then to wait for future SNe to occur. We observed approx- 
imately 160 galaxies with typical exposure times of 460s in 
F450W, F814W and F606W. These set of images produced 
pre-explosion environments for SNe 2003jg, 2004A, 2004dg, 
2005V, 2006ov, 2007aa, 2007gr and SN2008ax. The progen- 
itors of S N2004A and SN2008ax have been dete cted (see 
Sect l5T0l iHendrv et aT]|2006l ; ICrockett et~ail 120081 ) and the 
limits on the non-detections have proved very useful. This 
provides a benchmark for any future studies dedicated to 
detecting progenitors retrospectively. The depth reached in 
these images (typically mF8i4w — 25, 3<r) at a distance of 
~20 Mpc results in sensitivities down to Mfsi4w ~ —6.8 
(for a typical E[B -V)~ 0.2). 

The most fruitful method of detecting progenitors in 
the future would be to carry out a deep, wide-field survey 
of starforming galaxies within about 10 Mpc with the re- 
vived ACS and the new WF3 on HST after servicing mis- 
sion SM4. Reaching AB magnitudes of around 26 would be 
required to ensure the images go deep enough to detect pro- 
genitors down to around 8M©. For example SN2003gd was 
discovered at V = 25.8 (My = -4.4) and SN2003gd, 2005cs 
and 2008bk were all discovered at Mi ~ -6.5 ± 0.2 (all 
Vega based magnitudes). By restricting the distance limit 
to 10 Mpc would of course result in a lower rate of CCSNe 
(1-2 per year) in the sample of a round 100 major starforming 
galaxies (|Kennicutt et al]|200Sl ). but the discovery potential 
would be excellent. One further issue is the size of the star- 
forming disks of these galaxies which are significantly larger 
than the 3.4 arcmin diameter of ACS Wide-Field-Channel. 
The limited size of the WFPC2 and ACS cameras has re- 
sulted in many galaxies having been observed by HST before 
a CCSN occurs but the position falling outside the cam- 
era footprints (43 per cent : see Sect l2.2[) . Hence a carefully 
planned future survey of nearby galaxies must use multiple 
HST pointings to cover the full optical extents of the galaxy 
disks. A campaign of this extent is already in the arena of 
"multi-cycle treasury proposals" . 

A systematic survey beyond 10 Mpc that would signif- 
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icantly increase detection probability is probably not the 
best strategy. A distance limit increase to rioMpc (where 
)"iOMpc is the distance limit in units of 10 Mpc) would in- 
crease the number of galaxies by a factor of approximately 
100(ri OMpc — 1) and the time required would increase more 
than linearly (depending on the rioMpc distribution and 
galaxy sizes). The latter is not feasible in any reasonable 
allocation of time and given the difficulty in detecting pro- 
genitors beyond 10 Mpc is not an optimal strategy. 

We have found it essential that astrometry of the SNe 
employs images of the resolution of the HST pre-explosion 
images to avoid large astrometric uncertainties and am- 
biguous (or spurious) detection of progenitors. Ground- 
based (natural seeing) images are normally not adequate 
for providing geometric transformations between pre and 
post-explosion images to the 10-12 milli-arcsec level re- 
quired. Between 2001-2006 our group used HST to take 
post-explosion images at the same resolution (or higher, 
depending the camera pixel scales) as the pre-explosion 
HST images. This proved to be essential as there were sev- 
eral announcements (in IAU Circulars and Central Bureau 
Electronic Telegrams) of erroneous detections of progeni- 
tors using ground-based astrometry dVan Dvk et all l2003al ; 
iRichmond fc Modiazlteooll iLi et alfoooti ). Using unchecked 
and absolute astrometry to identify progenitors on ground- 
based images has also lead to erroneous clai ms of progen- 
itor d etections (SN 2002ap and SN 2004dj: ISmartt et all 
120021 ; ILi et al.l 12004 IWeis et al.l |2004 ) Images with adap- 
tive optics systems on large ground based telescopes (e.g. 
the VLT, Gemini and Keck) routinely deliver diffraction 
limited images in the K— band. These images are of suffi- 
cient resolution to use instead of HST follo w-up images (e.g. 
iGal-Yam et aLlbOOTtlCrockett et alj|2008h Rapid analysis of 
images after SNe discoveries and reporting of possible pro- 
genitors before extensive analysis has been completed is un- 
derstandable in a competitive field and helps the community 
to prioritise potentially interesting events. However claims 
of spatial coincidence require detailed differential astrome- 
try between post and pre-explosion images to precisions of 
around 10 milli-arcsec and comprehensive error analyses in 
order to be acceptable. Our experience has shown if this rule 
of thumb is not followed then erroneous results often follow. 
The flexibility and rapid reaction time of ground-based AO 
systems make them an excellent facility to provide precision 
differential astrometry (at the 10 milli-arcsec level) , and this 
has been our preferred strategy since 2007. 

Finally, we are now at the point at which, for some pro- 
genitors, we can return to see if the st ars have disappeared . 
This has been done for only SN 1987A (|Gilmozzi et alj|l987h 
and deep late-time images of SN 2004et alerted us to likely 
spurious detection of the progenitor. This area of late-time 
observations may lead to further surprises, but it is essen- 
tial it is pursued to verify the results claimed here, and else- 
where, for 2003gd, 2005cs and 2008bk. 



10 CONCLUSIONS 

This paper presents a consistent and homogeneous analysis 
of the constraints on progenitor stars of II-P SNe, within a 
volume and time limited search. The work builds and en- 
hances the discoveries and limits presented in the literature. 



There are now enough data that a statistical study of pro- 
genitor properties is feasible. This represents the culmina- 
tion of a a 10.5 year search for progenitor stars and the 
conclusions are as follows: 

• We have compiled all SNe discovered within a strict 
volume and time limit and reviewed the SN types reported 
in the literature. This gives an very good estimate of the 
relative rate of CCNSe in the Local Universe, at metallicities 
between LMC and solar. 

• Of the 55 type II-P SNe found, 20 have HST pre- 
explosion (or excellent ground-based) images available to 
search for progenitor stars. We summarise the data pre- 
sented in the literature to date on these SNe and carry out 
a consistent and homogeneous analysis of all events. Three 
groups of events are discussed - those with probable single 
star progenitors detected, those with upper limits to their 
luminosities and masses, and those falling on compact, un- 
resolved coeval clusters 

• The masses and mass limits of each are determined us- 
ing our STARS evolutionary models and a statistical anal- 
yses presented of the final masses. A maximum likelihood 
analysis suggests that the minimum mass for a type II-P 
SN is m m i n — 8.5_j 5 Mq. This is consistent with current 
estimates of the maximum mass that will produce a white 
dwarf in young clusters. 

• We have not detected any progenitors above 16Mq . As- 
suming a Salpeter IMF, the most likely maximum mass for a 
II-P SN is rtimax = 16.5±1.5Mq .This is not particularly sen- 
sitive to the IMF slope within the typical variations known 
in the Local Universes of ±0.7. 

• We suggest that there is a discrepancy between this 
maximum mass and our knowledge of massive star evolu- 
tion. Red supergiants between 17 - 30M© are not detected 
as progenitors but are predicted by theory to exist and this 
is supported by stellar population studies. We term this dis- 
crepancy the "red supergiant problem" . It is unlikely to be 
due to IMF variations and possible explanations include the 
possibility that we have systematically under estimated the 
stellar luminosities and masses due to foreground extinction 
or that the gap is filled with the other flavours of type II SNe 
(e.g. II-L, Iln and lib). However neither of these solutions 
are supported by current data. We suggest that these ob- 
jects may be forming black- holes with faint, or non-existent 
SN explosions. 

• Although low-luminosity SNe are already known (e.g. 
1999br, 2005cs) we suggest that these events are not the 
missing SNe - rather our analysis supports the interpretation 
that they are stars at the minimum mass limit for SNe. 

• We review the information on extinctions toward these 
SNe and their progenitors and compare it with that toward 
red supergiant populations in the LMC. We suggest that 
high extinction toward the SNe progenitors is unlikely to be 
the cause of the lack of detections of massive supergiants. 

• The search for progenitor stars should continue for ev- 
ery nearby SN that has deep, high resolution pre-explosion 
images. In particular the missing high mass red supergiants 
should be sought in both optical, near infra-red and mid 
infra-red images. If the limit for black hole formation is as 
low as 17M© then it bodes well for surveys f or disappearing 
stars as suggested by iKochanek et al.l (|2008t >. 
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APPENDIX A: COMPLETE LIST OF SNE 
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PROGENITORS 
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Table Al: SNe discovered between 1998-2008.5 in galaxies with recessional 
velocities less than 2000 kms -1 . The HST FOV column notes if the host 
galaxy has been observed by HST prior to explosion and if the position of the 
SN is "in" or "out" of the camera field-of-view. 



Supernova 


Galaxy 




Type 


HST FOV 


Comments 


Core-collapse 












1998A 


IC2627 


1975.7 


IIpcc 




Pastorello et al. (2005) SN1987A-like 


1998S 


NGC3877 


1114.8 


Iln 




Fassia et al. (20011 


1998bm 


IC2458 


1809.7 


11 




Li et al. (1998). unknown subtype 


1998bv 


PGC2302994 


1793.7 


II-P 




Kniazev et al. (1998) 


1998dl 


NGC1084 


1298.2 


II-P 




Filippenko fc De Breuck (1998a) 


1998dn 


NGC 337A 


997.1 


11 




Gomez & Lopez (2000), unknown subtype 


1999B 


UGC7189 


1962.1 


11 




Filippenko et al. (T999), unknown subtype 


1999an 


IC 755 


1607.2 


II 


in 


Wei et al. (1999). unknown subtype 


1999bg 


IC758 


1530 


II-P 




Avani &: Yamaoka (1999), amateur LC 


1999br 


NGC4900 


1012.6 


II-P 


in 


Pastorello et al. (2004) Low luminosity II-P 


1999bw 


NGC3198 


850.9 


LBV 


out 


Filippenko et al. (1999) 9 


1999el 


NGC6951 


1704.2 


Iln 


out 


Di Carlo et al. (2002) 


1999em 


NGC1637 


615.2 


II-P 


out 


Leonard et al. (2002) 


1999eu 


NGC1097 


1065.8 


II-P 


out 


Pastorello et al. (2004) Low luminosity II-P 


1999ev 


NGC4274 


1089.5 


II-P 


in 


Garnavich et al. (1999), amateur LC. 


1999ga 


NGC2442 


1149.6 


II-L 




Peculiar II-L (Pastorello et al. in prep) 


1999gi 


NGC3184 


765.8 


II-P 


in 


Leonard et al. (2002) 


1999gn 


NGC4303 


1616 


II-P 




Pastorello et al. (2004) possible low luminositv II-P 


1999gq 


NGC4523 


360.3 


II-P 




Jha et al. (1999), amateur LC 


2000ch 


NGC3432 


778.6 


LBV 


out 


Wasmer et al. (2004) 


2000db 


NGC3949 


1022.8 


II-P 


out 


Pastorello et al. (2000), amateur LC 


2000ds 


NGC2768 


1620.5 


lb 


in 


Filippenko & Chornock (2000) 


2000ew 


NGC3810 


1061.3 


Ic 


in 


Gerardy et al. (2002) 


2001B 


IC 391 


1815.9 


lb 


in 


Chornock & Filippenko (2001): Tsvetkov (2006) 


2001X 


NGC5921 


1581.6 


II-P 




Tsvetkov (2006) 


2001ac 


NGC3504 


1668.9 


LBV 


out 


Matheson & Calkins (2001) 


2001ci 


NGC3079 


1340.5 


Ic 


in 


Filippenko & Chornock (2001) 


2001du 


NGC1365 


1407.8 


II-P 


in 


Smartt et al. (2003); Van Dvk et al. (2003b) 


2001fv 


NGC3512 


1505.5 


II-P 




Matheson et al. (2001). amateur LC 


2001fz 


NGC2280 


1720.9 


II-P 




Matheson et al. (2001), amateur LC 


2001gd 


NGC5033 


1080.6 


lib 


out 


Stockdale et al. (2007) 


2001ig 


NGC7424 


754.4 


lib 


out 


Rvder et al. (2006); Maund et al. (2007) 


2002E 


NGC4129 


1152 


11 




Matheson & Berlind (2002), unknown subtype 


2002ao 


UGC 9299 


1614.9 


Ic 




Pastorello et al. (2008): narrow He lines 10 


2002ap 


NGC 628 


685.7 


Ic 


out 


Mazzali et al. (2002) 


2002bu 


NGC4242 


736.6 


Iln 


out 


Ayani et al. (2002), LC appears like II-L 


2002hc 


NGC2559 


1371 


II-L 




Wood-Vasev et al. (2002). similar to 1979C 


2002hh 


NGC6946 


318.5 


II-P 


out 


Pozzo et al. (2006) 


2002ji 


NGC3655 


1525.4 


Ib/c 




Ruiz-Laouente et al. (2002) 


2002jz 


UGC 2984 


1529 


Ic 




Hamuy et al. (2002), somewhat uncertain 


2002kg 


NGC2403 


370 


LBV 


out 


Maund et al. (2006) 


2003B 


NGC1097 


1065.8 


II-P 




M. Hamuy (priv. communication) 


2003J 


NGC4157 


1002.7 


II-P 




Avani et al. (2003), amateur LC 


2003Z 


NGC2742 


1517.7 


II-P 




Utrobin et al. (2007), Spiro et al. in prep.; low luminositv II-P 


2003bg 


M-05-10-15 


1151 


Ic 




Soderbere et al. (2006). evidence of H during evolution 


2003bk 


NGC4316 


1325.5 


11 




Phillips et al. (2003), subtype uknown 


2003ed 


NGC5303A 


1841 


lib 




Leonard et al. (2003) 


2003gd 


NGC 628 


685.7 


II-P 


in 


Hendrv et al. (20051 


2003gm 


NGC5334 


1433.8 


LBV 


in 


Maund et al. (2006) 


2003hn 


NGC1448 


911.3 


II-P 


out 


M. Hamuy (priv. communication) 


2003ie 


NGC4051 


917 


II-P 




See Sect. 5.9. Harutvunvan et al. (2008) 


2003jg 


NGC2997 


914.1 


Ib/c 


in 


Howell (2003) 


2004A 


NGC6207 


1090.3 


II-P 


in 


Hendrv et al. (2006) 


2004C 


NGC3683 


1944.3 


Ic 




Matheson et al. (2004a) 


2004G 


NGC5668 


1672.7 


II 


out 


Elias-Rosa et al. (2004), subtype unknown 



continued on next page 



9 http:/ /etacar. umn.edu/etainfo/related/ 

10 SNe 2002ao and 2006jc have been termed Ibn as they show narrow He lines due to circumstellar He rich shells. We discuss in Sect. 12. ll 
why we use Ic rather than Ibn 
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Table IA1I continued 



Supernova 


Galaxy 




Type 


HST flag 


Comments 


2004am 


NGC3034 


487.1 


II-P 


in 


See Sect. 5.11 


2004ao 


UGC10862 


1811.6 


lb 




Matheson et al. f 2004b) 


2004bm 


NGC3437 


1379.6 


Ic 




Folev et al. (20041 


2004cm 


NGC5486 


1648.2 


II-P 




Connolly (2004), probable low-luminosity II-P 


2004cz 


E407-G09 


1399.8 


II-P 




Folev et al. (2004) 


2004dg 


NGC5806 


1439.1 


II-P 


in 


Elias et al. (2004), and Sect. 5.12 


2004dj 


NGC2403 


370 


II-P 


out 


Vinko et al. (2006): Zhane et al. (2006) 


2004dk 


NGC6118 


1645.4 


lb 




Filippenko et al. (2004) 


2004ep 


IC2152 


1687.1 


11 




Filippcnko et al. (2004), subtype unknown 


2004et 


NGC6946 


318.5 


II-P 




Misra et al. (2007); Sahu et al. (2006) 


2004ez 


NGC3430 


1731.3 


II-P 




Filippcnko et al. (2004), amateur LC 


2004fc 


NGC 701 


1727.2 


II-P 




Salvo et al. (2004) 


2004gk 


IC3311 


-79 


Ic 




Ouimbv et al. (2004) 


2004gn 


NGC4527 


1780.3 


Ic 


out 


1990B-like, SNIFS spectrum 


2004gq 


NGC1832 


1779.7 


Ib/c 


_ 


Modiaz & Falco (2005); Gal- Yam et al. (2007) 


2004gt 


NGC4038 


1570.4 


Ib/c 


in 


Gal- Yam et al. (2005); Maund et al. (2005) 


2005V 


NGC2146 


1158.4 


Ib/c 


in 


Taubcnbcreer et al. (2005) 


2005ad 


NGC941 


1519.6 


II-P 




M. Hamuy (priv. communication), probable II-P 


2005ae 


E209-G09 


860.9 


lib 


in 


FiliDDenko & Folev (2005) 


2005af 


NGC4945 


373.5 


II-P 




Pereyra et al. (2006) 


2005at 


NGC6744 


617.6 


Ic 


out 


Schmidt fc Salvo (2005) 


2005ay 


NGC3938 


1016.8 


II-P 




Gal- Yam et al. (2008) 


2005cs 


NGC5194 


702.1 


II-P 


in 


Pastorello et al. (2006), low-luminosity II-P 


2005cz 


NGC4589 


2283 


lb 


in 


Leonard (2005) 12 


2005kl 


NGC4369 


1248.2 


Ic 




Taubcnbcrgcr et al. (2005) 


2006bc 


NGC 2397 


1059.7 


II-P 


in 


Patat et al. (2006), amateur LC 


2006bp 


NGC 3953 


1285.3 


II-P 




Dessart et al. (2008) 


2006fp 


UGC12182 


1807.5 


LBV 




Blondin et al. (2006) 


2006jc 


UGC 4904 


1818 


Ic 




Pastorello et al. (2007); Folev et al. (2007): narrow He 13 


2006my 


NGC 4651 


912.1 


II-P 


in 


Li et al. (2007) 


2006ov 


NGC 4303 


1616 


II-P 


in 


Li et al. (2007) 


2007C 


NGC4981 


1686 


lb 




Blondin & Calkins (2007b) 


2007Y 


NGC1187 


1213 


Ib/c 




Folatelli et al. (2007) 


2007aa 


NGC4030 


1475.7 


II-P 


in 


Folatelli et al. (2007). our own LC 


2007av 


NGC3279 


1435.5 


II-P 




Harutvunvan et al. (2007), amateur LC 


2007gr 


NGC1058 


633.8 


Ic 


in 


Valenti et al. (2008) 


2007it 


NGC5530 


1045.3 


II 




Contreras et al. (2007), unknown subtype 


2007oc 


NGC7418A 


1938.3 


II-P 




Olivares & Folatelli (2007), 1999em like 


2007od 


UGC12846 


1810.5 


II-P 




Blondin & Calkins (2007a), 1999em like 


2007sv 


UGC 5979 


1377.0 


LBV 




Harutvunvan et al. (2007). our own data 


2008M 


E121-G26 


1979 


II-P 




Folatelli et al. (2008) 


2008ax 


NGC4490 


797 


lib 


in 


Pastorello et al. (2008) 


2008bk 


NGC7793 


60 


II-P 


out 


Sect 5.2, Mattila et al. (2008) 


2008bo 


NGC6643 


1778 


lb 




Navasardvan et al. (2008) 


Type la 












1998aq 


NGC3982 


1388.4 


la 


in 


Branch et al. (2003) 


1998bn 


NGC4462 


1717.8 


la 




Patat & Maia (1998) 


1998bu 


NGC3368 


940.2 


la 


out 


Hernandez et al. (2000) 


1998dm 


PGC 005341 


1881.9 


la 




Filippenko & De Breuck (1998b) 


1999by 


NGC2841 


830.9 


la 




Garnavich et al. (2004) 


2000E 


NGC6951 


1704.2 


la 




Valentini et al. (2003) 


2001dp 


NGC3953 


1285.3 


la 


out 


Avani et al. (2001) 


2001el 


NGC1448 


911.3 


la 


out 


Krisciunas et al. (2003) 


2001fu 


M-03-23-11 


1555 


la 




Matheson & Berlind (2001) 


2002bo 


NGC3190 


1396.1 


la 


out 


Benetti et al. (2004) 


2002cv 


NGC3190 


1396.1 


la 


out 


Elias- Rosa et al. (2008) 



continued on next page 



11 http://www.supernovae.net/sn2004/sn2004gn.jpg 

12 Although V v ir outside limit, the TF/SBF distance from lTonrv et al I fcOOOl) puts it within our distance limit and it will be included 
in the stripped progenitor study of Crockett et al. in prep 

13 SNe 2002ao and 2006jc have been termed Ibn as they show narrow He lines due to circumstellar He rich shells. We discuss in Sect. 12. l1 
why we use Ic rather than Ibn 
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Table IA1I continued 



Supernova 


Galaxy 




Type 


HST flag 


Comments 


2002fk 


NGC1309 


1986.7 


la 




Ayani & Yamaoka (2002) 


2003cg 


NGC3169 


1236.7 


la 


in 


Elias-Rosa et al. (2006) 


2003gs 


NGC 936 


1274.3 


la 


in 


Gonzalez et al. (2003) 


2003hv 


NGC1201 


1490.3 


la 


in 


Gerardv et al. (2007) 


2003hx 


NGC2076 


1967.1 


la 




Misra et al. (2008) 


2003if 


NGC1302 


1504.8 


la 




Matheson et al. (2003) 


2004W 


NGC4649 


1197.8 


la 


in 


Moore et al. (2004) 


2004ab 


NGC5054 


1706.6 


la 


in 


Matheson et al. (2004) 


2004ea 


M-03-11-19 


1791.3 


la 




Gonzalez et al. (2004) 


2005cf 


M-01-39-03 


1976.8 


la 




Pastorello et al. (2007) 


2005cn 


NGC5061 


1946.6 


la 


out 


Modiaz et al. (2005) 


2005df 


NGC1559 


1004.4 


la 


in 


Gerardv et al. (2007) 


2005ke 


NGC1371 


1272.2 


la 




Immler et al. (2006) 


2006E 


NGC 5338 


895.4 


la 




Alderins et al. (2006) 


2006X 


NGC4321 


1679.1 


la 


in 


Wane et al. (2008) 


2006ce 


NGC 908 


1345.3 


la 




Blackman et al. (2006) 


2006dd 


NGC 1316 


1542.2 


la 


in 


Salvo et al. (2006) 


2006mq 


ESO 494-G26 


779.4 


la 




Prieto (2006) 


2006mr 


NGC 1316 


1542.2 


la 


in 


PhilliDs et al. (2006) 


2007af 


JNGG5584 


1704.7 


la 


in 


Simon et al. (2007) 


2007bm 


NGC3672 


1821.5 


la 




Navasardyaii ii et_al i (2007) 


2007gi 


NGC4036 


1629 


la 


in 


Harutvunvan et al. (2007) 


20071e 


NGC7721 


1972.1 


la 




FiliDDenko et al. (2007) 


2007on 


NGC1404 


1716.7 


la 


in 


Voss & Nelemans (2008) 


2007sa 


NGC3499 


1762.2 


la 




Aenoletto et al. (2007) 


2007sr 


NGC4038 


1570.4 


la 


in 


Umbriaco et al. (2007) 


Unclassified 












1998cf 


NGC3504 


1668.9 




in 




1999gs 


NGC4725 


1206 




out 




Uncertain 












M85 OT2006-1 


M85 


957 


SN or OT? 


in 


Nature uncertain, See Sect. 2.4 


2008S 


NGC6946 


318 


Iln or OT? 


out 


Nature uncertain, See Sect. 2.4 


NGC300 OT2008-1 


NGC300 


100 


SN or OT? 


in 


Nature uncertain, See Sect. 2.4 



